The scale of gravity and the cosmological constant within a landscape 
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It is possible that the scale of gravity, parameterized by the apparent Planck mass, may obtain 
different values within different universes in an encompassing multiverse. We investigate the range 
over which the Planck mass may scan while still satisfying anthropic constraints. The window for 
anthropically allowed values of the Planck mass may have important consequences for landscape 
predictions. For example, if the likelihood to observe some value of the Planck mass is weighted 
by the inflationary expansion factors of the universes that contain that value, then it appears 
extremely unlikely to observe the value of the Planck mass that is measured within our universe. 
This is another example of the runaway inflation problem discussed in recent literature. We also 
show that the window for the Planck mass significantly weakens the anthropic constraint on the 
cosmological constant when both are allowed to vary over a landscape. 



I. INTRODUCTION 

Theoretical results from inflationary cosmology [l| and 
from string theory 0, 0] motivate the possibility of an 
eternally inflating multiverse that is populated by an in- 
finite number of sub-universes, each obtained via local 
tunneling, diffusion, and/or classical slow-roll into one of 
a myriad of allowed meta-stable states. (For other mo- 
tivations to consider such a landscape see for example 
refs. [4(.) In this landscape picture, each of these uni- 
verses may contain different values for physical param- 
eters, or even different particles and interactions, than 
those that are observed within our local universe. How- 
ever, the anthropic principle 0] asserts that the physical 
laws that may be observed within any universe must be 
restricted to those that permit the evolution of observers 
in the first place. 

Combined with the anthropic principle, the landscape 
picture has emerged as a plausible explanation for many 
striking features of our universe. In particular, it has 
been used to justify the "unnatural" smallness of the 
Higgs mass [6j , the cosmological constant [7|, |8( , and the 
neutrino masses @ ; to predict the size of supersymmetry 
breaking 0, [l(| ; to describe the tilt in the spectrum of 
density perturbations [Tl| and other inflationary parame- 
ters [12| ; to constrain the baryon to pho ton ratio [13| and 
the ratio of baryons to dark matter [l4 , [To| , as well as to 
explain some seemingly fine-tuned relationships between 
parameters describing the theories of quantum electro- 
dynamics and quantum chromodynamics [161 ] . Neverthe- 
less, generating precise predictions from a landscape pic- 
ture faces several major challenges. 

One of these challenges is to identify an anthropic cri- 
teria that is both specific and compelling. Yet even af- 
ter such an anthropic condition has been defined, it is a 
daunting task to discern its environmental requirements, 
to deduce their implications for physical parameters, and 
then to derive the associated anthropic constraints. One 
may proceed by considering the variation in only one 
physical parameter, starting from its value within our 
universe. However, apparently tight constraints on any 



single parameter may be significantly weakened when 
more than one parameter is allowed to vary. This seems 
to be the case with both the Higgs mass 1171 fsee however 
ref. [l||) and the cosmological constant [l9l [2fj|. More- 
over, all anthropically allowed universes may not be con- 
nected by the continuous variation of physical parame- 
ters. For example, the seemingly viable "cold big bang" 
universe (2lj results from independently varying several 
cosmological parameters to values very far from those 
obtained within our universe. 

In addition, to calculate the expectation values of phys- 
ical parameters within a landscape requires determining 
an appropriate measure to weight among the various pos- 
sible universes [22l l23l . l24l [2o| . That is, a precise an- 
thropic criteria does not account for all of the selection 
effects that contribute to the probability for a particu- 
lar universe to be observed. Universes may be more or 
less likely based on how readily they are obtained via the 
physical dynamics that govern the multiverse. To ac- 
count for this requires a complete understanding of the 
multiverse, its landscape, and the governing theory. In- 
deed, there is a more subtle challenge underlying this 
program, which is to develop an appropriate and self- 
consistent calculus to regulate calculations involving the 
infinite number of infinitely expansive universes that may 
be contained within the multiverse [22|, |23|, |25| • 

Nevertheless, a set of hypotheses to resolve these chal- 
lenges may be excluded if it predicts a very low likeli- 
hood to observe a universe with some physical character- 
istic that our universe possesses. For example, consider 
a proposal that includes a specific notion of observer, a 
consistent calculus to determine expectation values over 
the multiverse, and a theory to describe the landscape 
including how meta-stable states are mapped onto uni- 
verses within the multiverse. If this proposal then pre- 
dicts that an exponentially small number of observers 
measure a cosmological constant at or below the value 
obtained in our universe, then the proposal and the spe- 
cific landscape in question are probably not both correct. 

We investigate the possibility that the scale of grav- 
ity may scan over the landscape. This is consistent with 
the results of ref. [26| . where only parameters with mass 
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dimension were found to vary over a model of the land- 
scape. We everywhere parameterize the scale of gravity 
using the (reduced) Planck mass mp. Although mp is 
commonly taken to be a fixed fundamental scale, this 
need not be the case. For example, the multiverse may 
be governed by a low-energy theory with the Lagrangian, 

C= l -^—gM 2 F[cj>]R + C^+C l , (1) 

where the fundamental mass scale is M, R is the Ricci 
scalar, is the collection of fields that specify the meta- 
stable state of a universe, is the effective Lagrangian 
for these fields, and £ is the effective Lagrangian for 
matter. We assume that within each meta-stable state 
the fields <f> are very massive and fixed to values <f> — > (/>* 
and are therefore non-dynamical. The strength of gravity 
will therefore be a constant within each meta-stable state, 
determined by the effective Planck mass, 



mp = y/F[4u]M. (2) 

We also assume that within each universe mp is fixed 
prior to slow-roll inflation. 

Alternatively, our analysis may be viewed in the so- 
called Einstein frame where the scale of gravity is ev- 
erywhere fixed. To accomplish this, one simply per- 
forms the conformal transformation g^ — > F^ 1 [(jj^g^v 
Then the scale of gravity is everywhere M, but all other 
parameters with dimension mass are scaled by the fac- 
tor F[0*] -1 / 2 . Thus our analysis is equivalent to fixing 
the Planck mass to be the fundamental scale of physics 
but varying all other mass scales uniformly. Stated an- 
other way, in our analysis the frame-independent ratio of 
masses m/mp scales as -Ff^*] -1 / 2 , where for example m 
may be the cutoff of the theory, the Higgs mass, or the 
scale of strong dynamics. This picture was pre viously 
suggested in the penultimate section of ref. [17]. The 
idea that the effective value of mp may vary across the 
multiverse within the context of Brans-Dicke theory was 
studied in refs. [IHtH^I- In addition, ref. [27j studied a 
model of the form of Eq. |1} to show how inflationary 
dynamics can explain the hierarchy between the appar- 
ent Planck scale and the electroweak scale. As described 
below, our focus is different from the focus of this work. 

We calculate the range over which mp may scan while 
still satisfying anthropic constraints. For completeness 
we consider a wide range of environmental constraints. 
These relate to halo, galaxy, and star formation, in addi- 
tion to galactic and stellar dynamics. We restrict atten- 
tion to universes that possess the same particles, inter- 
actions, couplings, and physical scales that are observed 
within our universe. Note that this means the cut-off 
to the low-energy effective theory depends on M, not 
nip. In addition, we assume that whatever mechanisms 
drive inflation and provide bareogenesis are unchanged 
(except insofar as they depend on mp) across the land- 
scape that we consider. Within the context of the model 



described by Eq. ([T]) , these assumptions require the exis- 
tence of a large number of states that have approximately 
equal particle physics parameters yet different values of 
mp. We do not explore the interesting case of a land- 
scape model that permits only correlated changes in mp 
and the other particle physics parameters. Note that the 
above restrictions are conservative in the sense that lift- 
ing them can only expand the range of allowed mp. 

Many anthropic constraints relate to the formation 
of galaxies and depend on the spectrum of energy den- 
sity fluctuations evaluated at matter-radiation equality. 
On any given distance scale, this spectrum has an ap- 
proximately Gaussian distribution about some root mean 
square (rms) amplitude, and we calculate constraints as 
if all fluctuations have the rms amplitude. Perhaps not 
surprisingly, the range of mp that is consistent with all 
of the anthropic constraints is rather narrow. It also de- 
pends on what models are chosen for inflation, baryoge- 
nesis, and the dark matter. As an example, if we assume 
that inflation is chaotic with potential V((p) = ^m 2 ip 2 , 
that baryogenesis results from efficient leptogenesis, and 
assume a weakly interacting massive particle (WIMP) to 
be the dark matter, then anthropic considerations com- 
bine to constrain mp to be 0.1 < mp < 1.5, where mp is 
the ratio between mp and the value obtained within our 
universe, about 2.4 x 10 18 GeV. 

Even a very narrow anthropic range for mp may have 
significant consequences for proposals to calculate its ex- 
pectation value within a landscape. In particular, it is 
plausible that some proposals will ultimately weight uni- 
verses in part according to their inflationary expansion 
factor. This expansion factor depends exponentially on 
the number of e-folds of inflation that the universe un- 
dergoes, which in turn depends on mp. In this case the 
probability distribution for mp will be peaked only where 
some other selection effect cancels this strong exponential 
dependence. This other selection effect could be a very 
sharp peak or boundary to the underlying landscape dis- 
tribution; otherwise the effect must come from an expo- 
nentially strong anthropic dependence on mp. Yet such 
a strong anthropic dependence on mp would be in con- 
flict with the observation that mp has even a narrow 
anthropic window in our universe. Thus we are forced 
to conclude that under these weighting schemes the ob- 
servation of our universe is either extremely atypical or 
our value of mp sits at some sharp peak or boundary 
in the underlying landscape distribution. This point is 
completely analogous to the "a-problem" and "Q catas- 
trophe" identified in refs. [28| . 

We here note that a runaway problem associated with 
varying the effective Planck mass during eternal inflation 
has already been discussed in refs. [23, 29]. These papers 
studied the evolution of mp in Brans-Dicke theory when 
the Brans-Dicke field is allowed to be dynamical during 
inflation. On the other hand, we study the case where 
the fields <f> in Eq. ([T]) are very massive and therefore 
non-dynamical. The difference between these scenarios 
is subtle because at some level <j> must be dynamical in 
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order for the landscape to be populated within the mul- 
tiverse. Our approach is to assume that the fields <f> are 
only dynamical at the very high energies that dominate 
the dynamics of the multiverse. At these energies the dy- 
namics of <j> could be described as in refs. [H,l2§| or they 
could be described by different effects. We simply treat 
these dynamics as unknown except to assume that the <p 
are fixed prior to the slow-roll inflation that eventually re- 
heats into each of the anthropically favorable low-energy 
universes such as our own. 

We also consider the anthropic window for the cosmo- 
logical constant A when both A and mp are allowed to 
(independently) scan over the landscape. Even when the 
allowed range for mp is relatively narrow, it still allows 
for a significant broadening of the allowed range for A. 
To see this, note that A is constrained only by Weinberg's 
anthropic bound Q, 



PA < Peq^e 



(3) 



Here p\ is the energy density in cosmological constant, 
p oq is the matter density at matter-radiation equality, 
and tr e q is the typical fluctuation in matter density at 
equality. The broadening occurs because for WIMP dark 
matter, decreasing mp significantly increases p oq and for 
most models also significantly increases a cq . For exam- 
ple, if we again assume chaotic inflation with potential 
V(ip) — ^m^ip 2 and that baryogenesis results from effi- 
cient leptogenesis, then pa may be over a million times 
the value observed within our universe when mp > 0.1. 
Of course, a larger anthropic window for p\ does not 
necessarily imply that our value of pA is less likely to be 
observed. We illustrate the distribution of observed val- 
ues of pa with a very simplified calculation. The results 
of this calculation suggest that to observe the cosmolog- 
ical constant at or below the level obtained within our 
universe is very unlikely unless the landscape distribu- 
tion of mp is dominated by values very near to or larger 
than the value obtained within our universe. 

The remainder of this paper is organized as follows. In 
Section |TT] we calculate the range of mp allowed by an- 
thropic constraints in universes otherwise like ours. Con- 
straints come from a variety of cosmological processes 
and we summarize our results in Section Iff HI Then in 
Section UIT1 we argue that the value of mp that we observe 
is extremely unlikely if universes within the landscape are 
weighted by their inflationary expansion factor. In this 
section we also discuss some caveats to this argument. 
The scenario where both mp and A may vary across the 
landscape is discussed in Section IIVI Finally, we draw 
our conclusions in Section fVl 



II. ANTHROPIC CONSTRAINTS ON THE 
SCALE OF GRAVITY 

It is straightforward to organize the immediate effects 
of changing the scale of gravity when all other mass scales 
and couplings are kept fixed (this implies that the cut-off 



of the theory is also fixed). Then scanning the Planck 
mass corresponds to changing the proportionality con- 
stant between the Einstein and the stress-energy tensors, 



G — mp 



(4) 



In a homogeneous universe this simply changes the re- 
lationship between the Hubble rate H, its time rate of 
change H, and the energy (p) and pressure (P) densities, 
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This bears upon anthropic conditions because the Hub- 
ble rate determines when particle interactions freeze out 
of equilibrium. This affects the relative densities of, for 
example, matter to radiation and protons to neutrons. 

Of course the universe is only approximately homo- 
geneous. According to the present understanding, inho- 
mogeneities are generated by quantum fluctuations in at 
least one scalar field as it exits the Hubble radius dur- 
ing (nearly) de Sitter expansion in the early universe. If 
this is the case, then the Hubble rate also effects the size 
of the initial inhomogeneities. At late times, these in- 
homogeneities re-enter the Hubble radius and the scale 
of gravity takes on a new role. Then gravity provides a 
self-interaction to over-densities that may cause them to 
grow. Over-densities that grow too large become gravita- 
tionally bound and separate from the cosmic expansion. 
Within these structures, the expansion of the universe is 
inconsequential but the scale of gravity still determines 
the internal dynamics. 

We analyze the anthropic significance of these effects 
in chronological order, beginning with the effects on infla- 
tion. We then discuss baryogenesis, big bang nucleosyn- 
thesis, matter domination, structure formation, stellar 
dynamics, and finally the stability of stellar systems. The 
anthropic constraints are displayed in Fig. [T] and summa- 
rized in Section III HI The complexity of this analysis, 
along with the many uncertainties in our understanding 
of various cosmological processes, make a precise deter- 
mination of anthropic constraints impractical. Therefore 
we strive for approximations that capture the key effects 
of scanning mp. Usually, we determine the dominant 
scaling behavior of a quantity with mp and cite a pre- 
cisely determined value from our universe to determine 
the value in another universe. Unless otherwise stated, 
the values for cosmological parameters within our uni- 
verse are taken from the tables in ref. [H[ (note how- 
ever that we work in terms of the reduced Planck mass). 
Throughout this paper we denote the ratio of a quantity 
to the value that it obtains within our universe using a 
hat, for example 



mp 



mp 



2.4 x 10 18 GeV 



(6) 



Finally, we use units where Ti = c = ks = 1. 
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A. Inflation 

An early period of inflation is believed to have homog- 
enized our universe and yet provided the seeds of cosmic 
structure through the generation of small density per- 
turbations (for reviews of inflation see for example [30(). 
We parameterize these effects using the total number of 
e-folds of inflation TV and the Bardeen curvature pertur- 
bation (. In principle, both iV and ( are constrained 
by anthropic considerations. Meanwhile, for inflation to 
occur in the first place requires that the Hubble radius ex- 
pand at a rate slower than the speed of light. This effect 
is parameterized by requiring that the first slow-roll pa- 
rameter, ei = —H/H 2 , is smaller than unity. Finally, at 
some point inflation must end and the universe must re- 
heat to establish the initial conditions for the subsequent 
big bang evolution. We discuss the mp dependence of 
each of these below. 



1. Satisfying Slow-Roll for N e-folds of Inflation 

For inflation to occur in the first place requires that 
at some time ej < 1. When inflation is driven by the 
potential energy of a canonical scalar field tp, the first 
slow-roll parameter can be written 



Mechanism to generate ( e . 



mp dependence of ( e . 



nil 



(7) 



where V is the inflaton potential and the subscript on 
V v denotes differentiation with respect to ip. Although 
at first glance ej appears to increase with increasing mp, 
this can usually be compensated for by starting the in- 
flaton ip further up the potential. This is the case with 
each of the canonical inflationary models presented be- 
low. Therefore we assume that the occurrence of inflation 
in the first place does not significantly constrain mp. 

Meanwhile, the total duration of inflation is con- 
strained by the need to homogenize a universe large 
enough to allow for the formation of structure. This 
constraint, however, is very weak, since our observable 
universe appears immensely larger than is necessary to 
form a galaxy. Moreover, inflationary scenarios that pre- 
dict the observed value of ( typically allow for far more 
c-folds of inflation than are necessary to enclose our uni- 
verse. Therefore we assume that the anthropic bound on 
N does not significantly constrain mp. 



2. The Curvature Perturbation ( 

Anthropic constraints on the amplitude of ( stem from 
primordial black hole production, structure formation, 
and the stability of stellar systems and are described in 
Sections III Dl III El and lll Gl Presently we discuss the de- 
pendence ((mp) for future reference. For inflation driven 



Inflation with V(ip) oc <p p 
Natural inflation 
Hybrid inflation 
Ghost inflation 
End of inflation scenario 
Curvaton scenario 
Inhomogeneous reheating 



372=2 



-5/2 



Xm p 

X" 1 ™.- 1 or xmp' 



TABLE I: The dependence of ( e q on mp for a variety of mech- 
anisms to generate the curvature perturbation. The result for 
inhomogeneous reheating depends on the relative size between 
X and other mass scales in the Lagrangian, and can interpo- 
late between the two dependences given above. 



by a canonical scalar tp, the curvature perturbation on a 
co-moving scale with wave- vector k is 



((fc)cx 



1/3/2 



mpV<p 



(8) 



k=aH 



Vip are evaluated when the scale k exits 



where V and 

the Hubble radius. Anthropic constraints on ( apply to 
scales k < fc cq , where k cq is the wave- vector of the Hub- 
ble radius at matter-radiation equality. The potential V 
evaluated when these scales first exit the Hubble radius 
may depend on mp even when V(tp) does not. 

In our universe, ( does not change appreciably with 
fc, and we assume this holds at least approximately in 
other universes. Therefore we take ( w ((fc e q) = (eq over 
all scales of interest. The scale fc oq of matter-radiation 
equality itself depends on mp. However this dependence 
is logarithmic and its effect on ( cq is suppressed by the 
smallness of the slow-roll parameter, so we ignore it. Still, 
to solve for ( eq requires to choose a specific model for in- 
flation. Since there is no standard model of inflation, we 
must be content with a only a plausible range for the 
dependence on mp. We deduce this range by studying 
several of the most popular models of inflation. The re- 
sults for chaotic inflation [3l| with V(f) oc <p p , for hybrid 
inflation |32l|. for natural inflation [331 ] . and for ghost in- 
flation [3J| are listed in Table HI 

The curvature perturbation ( Qq may also be generated 
at the very end of inflation 35] or even much later as 
in the curvaton (36j and inhomogeneous reheating (37j 
scenarios. In each of these models a sub-dominant scalar 
X receives fluctuations while the fluctuations to the in- 
flaton are presumed to be negligible. The fluctuations 
in x are then transferred to radiation either at the end 
of inflation, during reheating, or during a phase transi- 
tion much later. In each of these cases, the amplitude of 
the curvature perturbation depends on the local vacuum 
expectation value (vev) of x- The dependence of ( cq on 
X and mp when each of these mechanisms operates ef- 
ficiently (that is, when any reheating occurs far out of 
equilibrium) is listed in Table [I] 

We note that these models are very flexible to an- 
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thropic selection. Although in principle x may be set 
by interactions such that it is fixed among the set of uni- 
verses we consider, in most cases % is a stochastic variable 
over these universes. When £ eq is generated at the end of 
inflation or via the curvaton or inhomogeneous reheating 
scenarios, then this implies that £cq is also a stochastic 
variable over different universes. Therefore in these cases 
there exist universes with far different mp but the same 
£ eq as in our universe, as well as universes with the same 
mp but different £ cq . 

For future convenience we write £ eq in the form, 

Cc q « mp" , (9) 

where according to our notation ( oq is the curvature per- 
turbation relative its value in our universe. The various 
models of inflation that we studied suggest that we should 
restrict a to the range < a < 3. However, the mecha- 
nisms that generate Ccq at the end of or well after inflation 
may generate a wide range of £ cq for a wide range of mp. 
Although it is technically possible that keeping inflation- 
ary parameters fixed but varying mp will cause the dom- 
inant contribution to the curvature perturbation to shift 
from one mechanism to another, this scenario should still 
be well approximated by the above guidelines so long as 
the variation in mp is not too large. 

3. Reheating 

The reheating of the universe after inflation is achieved 
by coupling the inflaton to other degrees of freedom. In a 
typical model, after inflation the inflaton rocks within its 
potential well and red-shifts like matter. The radiative 
decay products of the inflaton then dominate the energy 
density of the universe only after the Hubble rate falls 
below the decay width T/. At this time the energy den- 
sity in the inflaton is pi — 3mpr|. If Tj is independent 
of mp, then the reheat temperature scales like 

f RH = mj 2 . (fO) 

We have no empirical knowledge about reheating other 
than that Trh is above the temperature of big bang nu- 
cleosynthesis. However, if our universe is described by 
a grand unified theory, then Trh must be below the 
temperature of monopole production. In addition, if net 
baryon number is not generated during reheating, then 
Trh must be high enough to support the dominant mech- 
anism of baryogenesis. We comment on constraints like 
these in Section ITTT1 



B. Baryogenesis 

We parameterize the net baryon number of a universe 
with the ratio between the number density of baryons and 
the number density of photons: rj = n(,/n 7 . Although 



we see no direct anthropic constraints on the value of n, 
it will enter into the anthropic constraints described in 
Sections III El and III Gl Presently, we seek to parameter- 
ize the dependence of rj on mp for future reference. As 
with inflation, there is no standard model of baryogene- 
sis. Therefore we must again content ourselves with only 
a range for the dependence on mp, based on the most 
plausible mechanisms. For a summary of these see for 
example the reviews of refs. [38l |. 

Perhaps the most plausible mechanism to produce net 
baryon number is leptogenesis [39| . For example, net 
lepton number is rather easily obtained by the out-of- 
cquilibrium decay of a right-handed neutrino (RHN) . The 
resulting lepton asymmetry can then be converted into 
net baryon number by sphaleron transitions within the 
Standard Model [4(j. The value of 77 that results from 
leptogenesis depends on how far out of equilibrium the 
RHN decays. If RHN decay occurs far out of equilibrium 
then the resulting baryon asymmetry 77 is independent of 
mp. Otherwise, it scales roughly according to r\ oc nip . 
Note that baryogenesis via leptogenesis requires that the 
RHN, for example, be produced in the first place. We 
comment on this requirement in Section IIIII 

The Standard Model of particle physics itself gener- 
ates appropriate conditions for baryogenesis, when it is 
augmented by relatively light supersymmetric (SUSY) 
scalars to strengthen the electroweak phase transition. 
The process of electroweak baryogenesis is complex; yet 
interestingly it operates independently the scale of grav- 
ity. That is, although universal expansion is necessary to 
decrease the temperature of the universe and thus spur 
the electroweak phase transition, this process is relatively 
independent the rate of temperature change. Therefore 
electroweak baryogenesis gives a baryon asymmetry that 
is independent of mp [38| . 

Finally, we look at Affleck Dine (AD) baryogenesis [4l| . 
This mechanism takes advantage of scalar fields that pos- 
sess baryon or lepton number, as would exist in SUSY 
or a grand unified theory (GUT). During inflation these 
fields may acquire large vevs, and then the influence of 
baryon non-conserving interactions on their subsequent 
evolution may generate significant baryon number. There 
are many models to implement AD baryogenesis and each 
may give a different dependence on mp. We simply giv e 
the result for a set of scenarios described in ref. [42j |. 
where the scalar fields overlap a SUSY flat direction that 
is lifted by non-renormalizable interactions and a nega- 
tive induced mass term during inflation. To fit the re- 
sulting baryon asymmetry to observation requires differ- 
ent parameterizations when different non-renormalizable 
interactions dominate; however in each case the depen- 
dence on mp is given by rj oc m p 3 ^ 2 

For future reference, it is convenient to write the 
baryon to photon ratio in the form 

77 ^p". (II) 

From the above discussion, we expect < (3 < 3/2, with 
perhaps the most plausible values being (3 = or 3/2. As 
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was the case with inflation, if more than one mechanism 
contributes to the net baryon number then we still expect 
the 77 that results to be well approximated by the above 
guidelines. 



C. Big Bang Nucleosynthesis 

The process of big bang nucleosynthesis (BBN) pop- 
ulates the universe with light elements. In particular, 
the mass fractions of hydrogen (X) and helium (Y) are 
important for anthropic considerations described in Sec- 
tions III El and III Fl Since we are only interested in X 
and Y , we may take a very simplified view of BBN. 
Specifically, we assume that BBN generates appreciable 
concentrations of only hydrogen and helium-4. This is 
clearly appropriate within our universe, where the other 
products of BBN account for only about 0.01% of the 
mass fraction of the universe. Although this fraction may 
change significantly for differing values of mp, it would 
take a very large variation in mp for this change to be- 
come significant next to X or Y . A basic description of 
BBN can be found in ref. [H[. 

When we approximate BBN to result in only hydrogen 
and helium-4, we require only the ratio of neutrons (n) 
to protons (p) to deduce X and Y . Specifically, 



1 — n/p 

If, Y 

1 + n/p 



1 - X. 



(12) 



A free neutron has AE w 1.3 MeV more energy than 
a free proton. Thus if neutrons and protons are kept 
in thermal equilibrium by interactions that convert each 
into the other, then n/p is given by the Boltzmann factor, 



n/p = exp(~AE/T) 



(13) 



where T is the temperature. Neutrons and protons are 
converted into each other via interactions such as n + 
v < — ► p+e, where the v denotes an electron neutrino and 
the e an electron. In fact, the rates of these interactions 
decrease faster with temperature than does the Hubble 
rate. Therefore, below some temperature Tp the universe 
expands too rapidly for, for example, an n and a v to 
find each other and convert into a p and an e. Below this 
temperature the relative concentrations of n and p are 
fixed and the interaction is said to freeze-out. 1 

The freeze-out temperature Tp is obtained by equat- 
ing the total rate of interactions converting neutrons to 



1 The decay of neutrons decreases n/p from its value at Tp by 
about 14% within our universe. This loss is determined by the 
time at which neutrons are efficiently captured into Helium, 
which depends on mp only via a logarithmic dependence on 
n(mp) l43l . Inspecting Fig. 4.4 of ref. [43ll . it can be shown 
that in the most extreme case of f} = 3/2, less than half of the 
neutrons decay for mp < 6. On the other hand, for mp > 6 
the helium fraction is less than a tenth of the hydrogen fraction. 
Therefore this effect is always negligible at our level of analysis. 



protons, r np , with the Hubble rate H. Then the neutron 
to proton ratio is n/p sw cxp (— AE/Tp). Our universe 
contains n/p « 1/7 such that X 3/4. Note that vary- 
ing mp so as to decrease Tp works to decrease n/p and 
therefore push X closer to unity. Since in our universe 
X sa 3/4, this effect is negligible at our level of analysis. 
On the other hand, for temperatures T > Tp the rate 
of conversion between neutrons and protons is r np cx T 5 
and is independent of n. Since H cx T 2 /mp, the freeze- 

— 1/3 

out temperature scales like Tp cx m P when Tp is larger 
than in our universe. Therefore we find, 



n/p: 



exp — ln(7) Ap 3 



(14) 



The fractions X and Y are given by Eqs. (fT2"|) . For ex- 
ample, when mp = 10, 5, 0.2, and 0.1 we have X w 1, 
0.9, 0.5, and 0.4, respectively. 



D. Matter Domination 

We have assumed that the early universe is dominated 
by relativistic degrees of freedom, at least since BBN. 
However as the universe cools, massive degrees of free- 
dom eventually become non-relativistic. As it becomes 
non-relativistic, the number density of this matter be- 
comes exponentially suppressed relative that of radiation. 
However, this dilution eventually causes matter to freeze 
out of equilibrium with the remaining radiation. Sub- 
sequently, the energy density of a massive species i will 
redshift as pi cx m^n.; while the energy density in radi- 
ation scales as p ra( j cx Tn ra d. Thus it is inevitable that 
matter should ultimately come to dominate the energy 
density of the universe (structure formation constraints 
ensure that the cosmological constant does not become 
significant before matter domination). 

Nevertheless, the energy density at matter-radiation 
equality and the fraction of matter in baryons are rele- 
vant to anthropic constraints described in Sections III El 
and III Gl In addition, our assumption that the early uni- 
verse is radiation dominated does not hold if £ is too 
large. In this case, primordial black holes may dominate 
the energy density of the universe while baryons are still 
relativistic. Then all of the baryons would be redshifted 
away or swallowed into black holes. This possibility is 
studied at the end of this section. In the following we 
neglect the neutrino content of the universe. Their influ- 
ence on cosmology is commonly viewed as insignificant 
and we do not expect this to change since as a hot relic 
their density relative baryons is fixed. In addition, we 
assume the dark matter to be a WIMP. This allows for 
relatively precise predictions, as opposed to, for example, 
axion dark matter where the density is set by a stochastic 
variable |44|. (Note however that the stochastic nature 
of axion dark matter makes this possibility more flexible 
to anthropic selection, see for example refs. [Til Il5|.) 

The energy density in a WIMP dark matter candi- 
date is set by the relic abundance that results from the 
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frcczc-out of annihilation interactions when the temper- 
ature drops below the mass of the WIMP. After matter- 
radiation equality this gives the scaling, 

p cdB1 oc nip 1 T 3 . (15) 

Meanwhile, after baryogenesis the relative abundance of 
baryons n is conserved. Therefore at temperatures below 
the nucleon mass the energy density in baryons scales like 

Pb oc nT 3 . (16) 

Using that in our universe Pb/pcdm ~ 1/5 and that the 
energy density in radiation scales as p ra d oc T 4 , we find 
the energy density at matter-radiation equality to be 

/5 q » Q?7 + j^p 1 ) • ( 17 ) 
Finally, we note the baryon fraction within matter, 

A * Q + fr 1 ™? 1 ) 1 ■ as) 

As described above, these results do not hold if £ 
is so large as to produce an abundance of primordial 
black holes (PBHs) 45] . Numerical analysis reveals that 
a PBH is formed when an ener gy d ensity fluctuation 
a > 0.7 enters the Hubble radius [46|]. Meanwhile, dur- 
ing radiation domination a = | £ at Hubble radius cross- 
ing [47l | . Therefore we require < 0.5 in order to pre- 
vent the formation of a PBH when the scale k enters the 
Hubble radius. In fact, this implies a somewhat stronger 
constraint on This is because £ is a stochastic vari- 
able with a Gaussian tail and because a PBH need not 
be formed each Hubble time in order for PBHs to domi- 
nate the energy density of the universe. This constraint 
is worked out in ref. [19| and we follow that analysis. 

The likelihood that a curvature perturbation with root 
mean square C is greater than or equal to 1/2 is /(£) = 
erfc(2 -3 / 2 £ -1 ), where the complementary error function 
is defined erfc(a;) = 2tt~ 1 / 2 e~ z dz. PBHs redshift 
like matter while the other degrees of freedom redshift 
like radiation. Therefore by matter-radiation equality 
PBHs will compose roughly (a oq /a p bh)/ of the energy 
density of the universe. Here a p bh is the scale factor at 
which newly formed PBHs have sufficient mass to persist 
until equality. Then PBHs do not dominate the energy 
density of the universe when 

erfc(2- 3 / 2 Ccq 1 )<^, (19) 

where we have neglected any tilt in £. This approxima- 
tion underestimates the largest £ e q since recent observa- 
tions suggest a negative tilt [48| . 

The erfc function depends very strongly on its argu- 
ment; therefore the nip dependence of the ratio a p bh/«eq 
is inconsequential for our analysis. As one looks to earlier 
times, the lifetimes of PBHs decreases more rapidly than 



cosmic time decreases. Therefore Eq. (fT9|) ensures that 
PBHs dominate at no time prior to equality. Solving 
for when PBH lifetimes equal about 70,000 years gives 
flpbh/flcq ~ 10~ 20 which gives Ccq < 6 x 10~ 2 . Translat- 
ing this into a constraint on nip, we find 

m£ > 7 x 10~ 4 , (20) 

where we have used that the density fluctuation at Hub- 
ble radius crossing is a » 5 x 10~ 5 [48]. Eq. ([20]) is always 
weaker than the stellar lifetime constraint of Section Ul Fl 

E. Structure Formation 

The formation of structure within our universe occurs 
in several stages. First, over-densities in the nearly pres- 
sureless dark matter begin to grow upon entry into the 
Hubble radius. During radiation domination, this growth 
is logarithmic with time, while after the dark matter 
comes to dominate the energy density of the universe 
over-densities grow in proportion to the growth in the 
cosmic scale factor. On the other hand, over-densities in 
the baryons cannot grow until after recombination. How- 
ever, within an e-fold or so after recombination they have 
grown to match the over-density in dark matter, and sub- 
sequently grow in proportion to the growth in the cosmic 
scale factor. When these over-densities have grown suffi- 
ciently they separate from the Hubble flow and virialize 
to form what are termed halos. 

After virialization, the cold dark matter within halos 
is stabilized against gravitational collapse by its inabil- 
ity to release its kinetic energy. However, the baryons 
within the halo must collapse beyond their initial virial- 
ization radius if they are to fragment and condense into 
galaxies and ultimately into stars. This requires that the 
baryons have a means to dissipate their thermal energy. 
The constraints on nip that are implied by these stages of 
structure formation are discussed in the sections below. 
Presently, we describe the initial growth in over-densities 
for future reference. The subject of galaxy formation and 
in particular star formation is complex and not yet fully 
understood. We rely heavily on the simplifying assump- 
tions and models of refs. [la. Ii~9j. 

We find it convenient to track the evolution of over- 
densities in position space, as opposed to Fourier space. 
At matter-radiation equality, the variance of energy den- 
sity fluctuations over scales with co-moving radius R is, 

r 00 dk 

<( R ) = V^ k2 W{kR)T(k)P{k) , (21) 
Jo Z7T 

where W is a window function that may be chosen to 
be a "top hat" with radius R, T is a transfer function 
to account for the evolution of perturbations between 
when they enter the Hubble radius and equality, and P 
is the primordial power spectrum of fluctuations, P(k) oc 
(£(fc) 2 ). We parameterize a co-moving scale with radius 
R according to the total mass /j, that is enclosed within a 
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sphere of radius R. In addition, we measure fi relative to 
the mass of our galaxy (more precisely the mass of our 
galaxy plus its dark matter halo); thus (j. = 1 corresponds 
to 10 12 M Q , where M Q is the mass of the sun. 
A numerical curve fit to Eq. (f2~Tj) gives [15j . 



1.45 x 10- d s(Ai)Coq(mp), 



(22) 



where the function s(/i) carries the scale dependence of 
tTeq- This scale dependence occurs because at the time 
of equality smaller scales have been within the Hubble 
radius for a longer time than larger scales. The function 
s(fi) is equivalent to Eq. (A13) in ref. [HI]. However, we 
have normalized s(p) such that s(l) = 1. In addition, 
we define the variable fi with respect to a different scale 
than the authors of ref. [la ] . Therefore within this paper 
s is given by, 



Km) = 



(0.76 ln[17 + ,u~ 1/3 ] - 0.22 



-0.27 



(23) 



Note that s is a decreasing function of ii. In addition, we 
emphasize that cr oq is the root mean square (rms) value of 
a Gaussian random field. Therefore constraints involving 
CT eq (or C eq ) are never sharp in the sense that they may 
be overcome by fluctuations that happen to be larger or 
smaller than is typical. 

After recombination but before the domination of cos- 
mological constant, a linear over-density is given by @ 



3 a 

5 a cc 



' cq • 



(24) 



Soon after recombination the first term is negligible. 
Eq. (|2"4|) is accurate until cosmological constant domi- 
nation, after which er grows by another factor of about 
1.44 and then stops. An over-density separates from the 
Hubble flow and virializes when a linear analysis gives 
a = 1.69 [4!|. Thus the cosmic mean energy density at 
virialization is, 



3 a. 



cq 



5 1.69 



Pcq 



1.4 X KT^eqCeq 



3 s 3 . 



(25) 



The energy density within the condensed halo is larger 
by roughly a factor of 187r 2 . We denote this as, 



18ttV ~ 2.4 x 10" 



3 D 3 



PeqCeq S 



(26) 



Note that these quantities depend on both mp and the 
mass scale \i of the virialized halo. 

The above description of halo formation relies on three 
important aspects of the standard cosmology: we assume 
that the dark matter density dominates over the baryon 
density, that recombination occurs before the virializa- 
tion of the dark matter halo, and that virialization oc- 
curs before the domination of cosmological constant. En- 
forcing the above conditions implies constraints on mp. 
These are discussed in the next section. We discuss the 
possibility for a non-standard path toward structure for- 
mation in Appendix [A"l 



1. Halo Virialization 

Before proceeding to galaxy formation, we must ensure 
that over-densities separate from the cosmic expansion 
and virialize before the domination of the cosmological 
constant halts their growth. As mentioned above, an 
over-density a has separated from the Hubble flow when 
a linear analysis gives a > 1.69. On the other hand, 
the maximum size that is reached by an rms linear over- 
density is given by, 



1.44 X^a, 
5 a c 



eq 



3.20 



(27) 



where a\ is the scale factor at p\ domination. Therefore 
the requirement that <Joo > 1.69 gives 



Pcq< 



Ccq* 3 >0.1. 



(28) 



Since with a larger value for mp structures form later, 
this is a constraint against increasing mp. Substituting 
previous results into Eq. (|28p gives 



6 Wr 



5 \ 4 
+ -m p 1 J mp 3a s(/i) 3 > 0.1 



(29) 



The curves that saturate this inequality for various 
choices of a and f3 are displayed in Fig. Q] under the label 
Voo = 1-69." In light of these plots, our value of mp 
may be construed as nearly saturating this constraint. 
However, this perception derives from the strong mp de- 
pendence of pvir, and not from our galaxy being at the 
edge of saturating the Weinberg bound. In addition this 
constraint, analogous to many others below, holds for an 
rms fluctuation a eq but is weaker for larger fluctuations. 



For these reasons we emphasize that although Eq. 
does not allow for mp to be increased significantly, it is 
still true that our value of mp is not at the edge of the 
anthropic range. Since the existence of an anthropically 
allowed window surrounding our value of mp is essential 
to the arguments of Section IIII1 we provide a more elab- 
orate discussion of this boundary in Appendix |5J This 
may serve as an illustration of how 'soft' are other con- 
straints that depend on cr eq . 

The above analysis assumes that halos virialize at least 
an e-fold or so after recombination. This is to ensure 
that baryons may collapse into the dark matter poten- 
tial wells and participate in the virialization. Thus we 
require p*/p rC c ^ e -3 . The energy density at recombina- 
tion is set by the temperature of recombination, which 
depends only logarithmically on mp and r\. We ignore 
this logarithmic dependence and take T rec ~ 3000 K in 
every universe that we consider. Using that at any time 
after equality oc p roc cx T 3 , we find for this constraint, 



Pcq^Ccq^^^SxlO 



-8 



(30) 



Since with a smaller value of mp the matter energy den- 
sity at equality and the amplitude of density perturba- 
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tions are both larger, this is a constraint against decreas- 
ing top . In terms of rhp and p, this gives 



6™* 



-3 



,3o 



s([i)- 3 > 8 x 1CT 8 . (31) 



The curves that saturate this constraint for various a and 
(3 are displayed in Fig. Q] with the label "i roc = e _2 / 3 t vir ." 

Finally, we require that dark matter dominate over 
baryonic matter so that dark matter potential wells are 
deep enough to condense baryon over-densities after re- 
combination. This simply translates into the constraint 
fb ;$ 1/2, which gives 



m P < b 1/{1 - f3) for < 1 
mp £ S 1 ^ 1- ^ for (3 > 1 . 



(32) 

For (3=1 this argument provides no constraint on mp 
since in that case fb is independent mp. 



2. Galaxy Formation 

Although the dark matter within a halo cannot dis- 
sipate its kinetic energy to further collapse, the baryons 
may do so via electromagnetic interactions. If the cooling 
timescale r coo i is less than the timescale of gravitational 
dynamics r grav , then not only do the baryons collapse, 
but perturbations in the baryon density fragment into 
smaller structures. These structures ultimately fall into 
a rotationally supported disk. Perturbations may further 
fragment if the disk satisfies the Jeans instability crite- 
ria, which is that T grav be less than the time it takes for 
a pressure wave to traverse the perturbation. Fragmen- 
tation continues until perturbations become Jeans-stable 
and over-densities relax adiabatically into hot balls of 
gas. This appears to be the path by which halos within 
our universe ultimately condense in to g alaxies of stars 
(for background see for example ref. [Hof). 

In order to ensure galactic dynamics similar to those 
within our universe, one might therefore first impose that 



for typical halos, r c , 



< 



We take the dynamical 



timescale of the halo to be the time it would take for a 



2 We note that more careful considerations involving galactic dy- 
namics may suggest a far weaker constraint than the one we 



pursue. We require r coo i < 



grav 



where each timescale is evalu- 



ated at virialization. However, the baryons within a halo will 
cool even if this condition is not satisfied. As described in 
ref. [5lll . this cooling pushes the gas of baryons along a curve 
in the temperature-density phase space that eventually leads to 
the condition t coo i < T grav being satisfied, albeit at a much later 
time. It is then necessary to consider an y o ther factors that 
might constraint the timescale t coo i. Ref. [15H has pointed out 
that if t coo i is not much smaller than the Hubble time, then 
baryons do not cool significantly before being reheated by halo 
mergers. Eventually halo merging ceases due to cosmological 
constant domination. However, even then, one must worry about 
too large a faction of baryons evaporating from the halo before 



test particle to free-fall to the center of the halo, 
spherical halo of constant density, this is 



'3 2 -1/2 



For a 



(33) 



For the cooling timescale we use the total thermal energy 
divided by the rate of energy loss, per unit volume: 



7"cool 



3 hp v 



2 m N p b A c 



(34) 



where pb is the mean molecular weight of the baryons 
in the halo (in units of the nucleon mass mjv), Ty^ is 
the mean temperature of the halo, and A c is the rate of 
energy loss per unit volume. The quantity fb Pvir/mNPb 
is the baryon number density, including electrons. 

The mean molecular weight depends on the ionization 
fraction and the hydrogen mass fraction of the halo. For 
example, for a fully ionized halo we have, 



Pb 



re H + 4n Hc ^ 4 

e + «H + "He 3 + 5X ' 



(35) 



where the subscripts denote electrons, hydrogen, or he- 
lium. Note that pb never strays more than a factor of two 
from unity. To estimate the temperature of the halo, we 
first note that in a virialized halo the mean kinetic energy 
equals half the mean gravitational binding energy. Thus 
for a halo of mass M we write Mv 2 ir w ■^m^, 2 M 2 R~ l , 
where V v \ r IS cl characteristic velocity for virialized parti- 

?3 



cles and R is the radius of the halo. Since M ~ 
^p b m N v 2 ir , 



and since T v \ r as ^UbiriMV 2 :,., we obtain 



^R A n 

3 Jt fvir 



1 

78 



PbrriNmp M '"p. 



2 S/f 2/3„l/3 
vir ' 



(36) 



Note that T v j r cx cr eq so that T v j r is a stochastic variable 
for halos of a given mass. 

The baryons within a halo may cool via Compton scat- 
tering, bremsstrahlung, the excitation of hydrogen or he- 
lium lines, in addition to other mechanisms. These all 
contribute to the rate of thermal energy dissipation A c . 
Thus A c is a complicated function of temperature, which 
also depends on the halo composition and therefore the 
hydrogen and helium fractions X and Y. Rather than at- 
tempt an estimate of A c , we use the cooling rates given in 
refs. These include the processes listed above, but 

we neglect the possibility for molecular cooling, which is 
insignificant at the temperatures we consider. The galac- 
tic cooling constraint t coo i ;$ T grav is now, 



m^AJ^Ar^^SxlO- 3 . 



(37) 



The curves that saturate this inequality are displayed in 
Fig. [1] with the label "r coo i = T gra v" 



they cool sufficiently to sink deeper into the gravitational poten- 
tial well [TH |. These considerations are beyond the scope of this 
paper. 
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As noted in ref. |l5|, if a galaxy contains too little 
mass then early supernovae may blow away a significant 
fraction of its baryons when they explode. We expect 
the effects of a supernova to be relatively localized if the 
gravitational binding energy of the galaxy by far exceeds 
the energy released in the supernova [15| . This can be 
ensured by requiring that the energy released in the su- 
pernova be less than the halo binding energy, 



M 2 



bind 



40ttto p R 



1 

26 



±m^M^p^ 



(38) 



Note that the baryons within a galaxy are much more 
tightly bound than when in the original halo (see for 
example the galactic disk estimates of Sections III E 31 
and III Gp . Therefore Eq. (f3"B")) is a significant underes- 
timate of the binding energy of a galaxy. 

We expect the energy released in a supernova to scale 
roughly as the binding energy of a Chandrasekhar mass 
at its Schwarzschild radius [151 ] , or as the binding energy 
of a typical star, both of which are proportional to mp 
(see Section Hi F|) . Thus we write this energy E sn mp, 
where E sn « 10 51 erg is the typical supernova energy 
within our universe. Requiring that £bind E sn gives, 



5 /Si /3 CcqM 5/3 s>4x 10- 



(39) 



Since halos of a given mass become more weakly gravi- 
tationally bound as mp is increased, Eq. (|3"9"|) is a con- 
straint against increasing mp. After inserting previous 
results this becomes, 



1 » -p , 5 » -l 
g TO p + g m p 



\4/3 

J m~ 5 ' a fi 5/3 s(fi) > 4 x 10~ 9 . (40) 



The curves that saturate Eq. 
under the label "-Ebind = E su 



40)) are displayed in Fig. [T] 



The mass of the disk is simply the mass of the baryons 
in the halo, M^sk = foM. Meanwhile, the peculiar ve- 
locity is related to the temperature of the disk. This 
temperature will be the lowest temperature to which the 
baryons can cool as they collapse, which is roughly set by 
the hydrogen line temperature Th « 10 4 K. 3 Therefore 
the typical peculiar velocity may be written, 



/ 3T H 



(42) 



The circular velocity v c is deduced by conserving angular 
momentum as the baryonic halo collapses. On the one 
hand, the disk angular momentum can be roughly written 
as -Rdisk-Mdisk^e- On the other hand, the baryons in a halo 
start with angular momentum 



^hXm^M^R^, 

/07T * 

where A is the dimensionless spin parameter [541 ] . 



A = V8n 



JEbi, 



mp 



(43) 



where J is the magnitude of the angular momentum, 
-Ebind is the gravitational binding energy, and all of the 
above quantities are evaluated for the original halo. Be- 
fore the gravitational collapse of the baryons out of the 
dark matter halo, it is reasonable to assume that the 
angular momentum of the baryons and dark matter are 
equally distributed according to mass, such that initially 
Jb = MbJh/Mh where the subscripts b and h refer to 
baryonic and dark matter halo quantities [53| . Then 
assuming that angular momentum is conserved as the 
baryons in the halo collapse, we find 



A 



(44) 



3. Star Formation 

If the above conditions are met, the baryons in a halo 
will radiate away energy and settle into a disk supported 
by its angular momentum. We then require that the 
disk fragment so that ultimately stars may form [l5| . 
The stability of galactic disks against both radial and 
vertical perturbations can be studied using a standard 
Jeans analysis, which compares the dynamical timescale 
Tgrav to the time it takes a pressure wave to traverse the 
perturbation. It turns out that the stability criteria for 
the two modes differ by only an order unity coefficient [H, 
IHol . [53l ] . In ref. [53| it is shown that for perturbations in 
the vertical direction, perturbations are unstable when 
the total mass of the disk satisfies, 



M disk > 120 mpv p v c R disk . 



(41) 



where v p is the typical peculiar velocity of particles in 
the disk, v c is the circular velocity of these particles, and 
i?disk is the disk radius. 



We do not require to solve for i?disk after this expression 
for v c is inserted into Eq. (|4"Tj) . 

Eq. (|4"Tj) can now be written in the simple form (l5j . 



A>4A 



T ■ 

J vir 



1/2 



(45) 



The spin parameter A is different for different halos, but 
it is roughly independent of mp, /1 and the amplitude 
of density perturbations, and it typically lies near A w 
0.05 Substituting into Eq. (g5D gives 



(46) 



^hm-'pl/^^s 1 / 2 > 0.2. 



3 Although the rate of galactic cooling is reduced below the tem- 
perature of hydrogen line freeze-out, cooling still proceeds via 
molecular transitions in, for example, H2- Therefore T m ; n may 
become very small if one is willin g to wait a long time before 
disk fragmentation. We follow ref. [15l | and study when galactic 
dynamics are similar to those within our universe. 
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This results in a constraint against increasing mp. We 
use previous results to convert this into a constraint on 
mp and //, which gives 



1 ~ -0 . u - -l 
-mp +-m P 



-1/3 



0.2, (47) 



where ^(mp) depends relatively weakly on mp and has 
been ignored. The curves that saturate this constraint 
for various a and /3 are displayed in Fig. [T] with the label 
"disk inst.." 

It is essential that at some point fragmentation ceases 
so that over-densities can smoothly collapse into a star. 
The process of fragmentation may be seen to terminate 
when individual fragments become sufficiently opaque so 
as to trap most of their radiation [551 ]. In the Jeans pic- 
ture, this allows the temperature of a perturbation to rise 
and correspondingly increase the sound speed and thus 
prevent further fragmentation. As described in ref. [55| . 
the mass scale at which this occurs is relatively indepen- 
dent the dominant contributions to the cooling rate and 
opacity and gives a typical stellar mass that scales like 
mp. Interestingly, this is the same scaling behavior that 
restricts the sizes of stars based on their internal tem- 
perature being high enough to fuse hydrogen and their 
radiation pressure being low enough so as to not blow 
the star apart. We elaborate on this in the next section. 



F. Stellar Dynamics 

We have so far ensured that the fragmentation of over- 
densities persists on all scales greater than a relatively 
small scale that is roughly proportional to mp. We now 
require that the temperature within some of the remain- 
ing structures is sufficient to fuse hydrogen to form a star. 
It is possible that the mere existence of stars is not a suf- 
ficient condition for the existence of observers. Therefore 
we also consider the requirement that some of these stars 
supernova in order to generate heavy elements. In ad- 
dition, we consider the requirement that some of these 
stars have both a surface temperature within a factor of 
two that of the sun and a main-sequence lifetime of at 
least a few billion years. Our motivation for selecting 
these specific criteria is simple. Without knowing what 
are the necessary conditions for observers to arise within 
a stellar system, we study what seem at least to be two 
sufficient conditions. 

Ref. [HI has studied the basic requirements that con- 
strain the properties of stars. For decreasing stellar mass, 
the central temperature must be above some minimum 
temperature T nuc that is necessary to fuse hydrogen. The 
central temperature within a low-mass star is estimated 
by balancing the influences of gravitational pressure and 
electron degeneracy. T his g ives for the central tempera- 
ture of a low-mass star [16| , 



where M is the stellar mass. The least massive stars 
have T c = T nuc . Since T nuc is independent of mp, in any 
universe these stars have mass 



M min cx m P 



(49) 



On the other hand, there is also an upper limit to the 
mass of a star. If the radiation pressure within a star 
well exceeds the gravitational pressure that its mass can 
provide, then the star itself becomes unstable upon the 
ignition of its core. This constrains the maximum mass 
that a star may have to satisfy the scaling [l6| , 



M max cx m P . 



(50) 



Note that the minimum and maximum mass of a star 
both scale as nip. This is also the scaling of the typi- 
cal mass that becomes sufficiently opaque to prevent fur- 
ther fragmentation. This means that, given a window of 
masses for which stars exist in our universe, there will 
also be such a window within universes with significantly 
different values of mp. 



1. Stellar Lifetimes and Spectra 

We now seek constraints to ensure that some of the 
stars produced within a particular universe have surface 
temperatures and main-sequence lifetimes that appear 
to be sufficient for the evolution of observers. Our pur- 
pose in investigating this condition is to ensure that we 
do not overlook what might be viewed as an important 
anthropic constraint. Therefore we adopt a very restric- 
tive perspective and require that some stars have surface 
temperatures of at least 3500 K and that these stars have 
main-sequence lifetimes greater than the timescale of bi- 
ological evolution, t cvo i w5x 10 9 yrs. This surface tem- 
perature is chosen in part to simplify the calculation of 
stellar lifetimes and in part because a black body at this 
temperature radiates a significant fraction of its power 
into the frequency band accessible to chemistry. The 
evolutionary timescale t cvo i should be understood to in- 
clude the time required for a planet to condense and cool, 
minus the time it takes the star to reach main-sequence 
hydrogen burning. This time may be different for differ- 
ent planets, but we do not expect it to form the dominant 
contribution to t cvo i. 

The main-sequence lifetime of a star is roughly equal 
to the available energy of the star divided by the typical 
rate that it radiates energy away, 



r* cx XML' 



(51) 



T c cx mp 4 M 4/3 . 



(48) 



Here L is the typical luminosity of the star during main- 
sequence and X is the hydrogen fraction. We assume that 
differences in composition other than differences in the 
hydrogen fraction have benign consequences. In addition, 
we neglect the mp dependence of the mean molecular 
weight /z*, since changes by only roughly a factor of 
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two as X ranges from zero to unity. For an introduction 
to the major concepts of stellar astrophysics, we have 
found useful refs. [Eg, [58[ . 

Our narrow purpose allows for a simplified analysis of 
the necessary stellar dynamics. Since we specify stars 
by their surface temperature, we write the luminosity 
L oc R 2 Tf for stellar radius R and surface temperature 
T s . To eliminate R, we note that the central temperature 
of a star scales as 



Thus we can write, 



T* OC 



1 M 



XmpT 2 



(52) 



(53) 



The lifetime of a star is maximized by considering the 
minimum allowed surface temperature, in this case T s = 
3500 K. Since both the lower end and the upper end of 
the window for stellar masses scale cLS 7Tip , clS 0i basic 
approximation we may take this to be the scaling for all 
stellar masses at fixed T s and T c . Combining this with 
Eqs. ((SI]), dSU), and ([53]) gives, 



Xm\ 



(54) 



where f* is measured in units of the main-sequence life- 
time of these stars within our universe. This lifetime is 
roughly 100 billion years [5{|. Therefore the constraint 



> 



becomes 



Xrhp > 5 x 10" 



(55) 



This requires that mp satisfy mp ^ 0.1. 

We now consider this analysis in a little more detail. 
In particular, we consider the effects of convection and 
electron degeneracy explicitly in order to motivate that 
we can keep T s /T c fixed while scaling mp and that the 
stellar mass scales like mp at fixed T c . An ionized star in 
which radiation pressure can be neglected and in which 
the energy transport is everywhere dominated by convec- 
tion is well approximated as a polytrope with polytropic 
index i — 3/2 [111 HI]. This implies certain scalings be- 
tween stellar properties and in particular that for these 
stars T s /T c is independent of mp, M and R. Within our 
universe, stars with surface temperatures at and below 
3500 K have masses M < 0.35M Q and are well approxi- 
mated by these polytropes. In addition, it can be shown 
that stars defined by these temperatures remain convec- 
tion dominated as mp is decreased [58[ . To see this intu- 
itively, note that convection is driven by tidal forces and, 
all else being equal, one expects the tidal forces within a 
star to increase as mp is decreased. Therefore, we expect 
T s /T c to be fixed for stars with T s — 3500 K as mp is 
decreased. 

Meanwhile, so long as electron degeneracy is significant 
within the center of the star, the mass required to achieve 
a fixed central temperature scales like mp [l6[. Stars 



with a surface temperature of 3500 K are indeed partially 
degenerate within our universe, but they could become 
non-degenerate after some amount of scaling M oc mp. 
To see that there can exist stars for which the degeneracy 
remains fixed, consider again the polytrope model. The 
electron degeneracy at the center of a star is a function 

of the ratio n e T c 3 ^ 2 , where n e is the number density of 
electrons. The scaling relations applicable to an i = 3/2 
polytrope imply that the electron density n e scales like 
the average density of the star, which at constant cen- 
tral temperature scales like m P M~ 2 . Next note that the 
electron degeneracy at the center of a star is a constant 
if the mass of the star scales as M oc mp for fixed cen- 
tral temperature T c . This is precisely the scaling that 
describes a partially degenerate star, which means that a 
partially degenerate star remains partially degenerate as 
mp is scaled while keeping T c fixed. Therefore we expect 
stars with a surface temperature of 3500 K to remain 
partially degenerate as we decrease mp keeping T s fixed, 
so that indeed M oc mp. These arguments justify the 
constraint Eq. Q55[) . 

It is illuminating to consider a different form of analy- 
sis. This applies to ionized stars where radiation pressure 
can be neglected but radiation dominates over convection 
in the transport of energy. Then the scaling of the stel- 
lar mass M with mp for fixed T c may be obtained for 
a class of stars (so-called "homologous stars" 4 ) by ap- 
plying homologous transformations to the equations of 
hydrodynamical equilibrium. To perform such an anal- 
ysis, we must phenomenologically model the opacity of 
the star and the rate of energy generation per unit mass 
with the respective formulae [561 ], 



e « e X 2 P T v . 



(56) 



The terms kq and eo are constants while the exponents 
n, s, and v depend on the physical properties of the star. 
The Kramers opacity of intermediate mass stars such as 
our sun is modeled using n = 1 and s = 7/2. Then it 
can be shown that, independent of v, the stellar mass 
scales like M oc X 1 / 3 mp ^ 3 for fixed central tempera- 
ture [56j]. For such stars T s /T c is not constant, so one 
cannot use Eq. ([53]). However, at fixed central tempera- 
ture the homology transformations determine the scaling 
of the luminosity to be I oc X 2 m P M" 1 oc X 5 / 3 mp 3 . 
Inserting these scalings into Eq. (|5"Tj) gives 



2/3 



(57) 



where is the main-sequence lifetime in units of an ap- 
propriate lifetime evaluated within our universe. 



4 This is a very restrictive class, since by definition the mass distri- 
bution for two homologous stars of mass Mi and radius Ri must 
satisfy m\(r)/M\ = m,2(r R2 / Ri) / 'M2 , where rrii(r) is the mass 
contained within a sphere of radius r. 
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This gives a slightly weaker dependence on mp than 
Eq. (|55[) ; however stars that are well-described by these 
approximations have shorter lifetimes than those de- 
scribed by an i — 3/2 polytrope. In addition, the physical 
characteristics that make these approximations applica- 
ble will not continue to describe stars as we decrease mp 
with T s fixed. (They do continue to describe stars as mp 
is increased for fixed T s , since this tends to lessen the im- 
portance of convection.) Nevertheless, this confirms the 
qualitative form of Eq. ([55]) and suggests that an analo- 
gous analysis would apply if surface temperatures closer 
to that of the sun were demanded. 



2. Heavy Element Production 

Supernovae are believed to be the exclusive source of 
heavy elements within our universe. However, the dy- 
namics of supernovae are very complex and are still not 
fully understood (for reviews see for example refs. [60l]). 
Therefore ensuring the existence of supernovae in uni- 
verses with differing values of mp is clearly speculative. 
In this section we simply provide some qualitative re- 
marks in support of this possibility. There are many 
types of supernovae within our universe; for convenience 
we focus on what are called type la supernovae. 

Type la supernovae are understood to erupt via the 
accretion of matter by a white dwarf star. Meanwhile, 
white dwarfs are created when a star has consumed all 
of its hydrogen and helium fuel but does not possess suf- 
ficient mass either to drive its central temperature high 
enough to ignite carbon fusion or to form a black hole. 
According to the scaling relationships discussed in the 
previous section, the first condition is always satisfied 
given that it is satisfied within our universe. On the 
other hand, the Schwarzschild radius and the physical 
radius of a star at fixed central temperature both scale 
as R cx rrip 2 M . Therefore white dwarfs will exist in all 
of the universes that we consider. 

The supernova of an accreting white dwarf proceeds 
when its growing mass reaches the Chandrasekhar limit 
and the star becomes unstable through the nuclear ig- 
nition of its carbon. The relevant physical scales for 
this phenomena are set by the Chandrasekhar mass and 
the binding energy of a Chandrasekhar mass at about a 
Schwarzschild radius. Since these and the typical stellar 
mass all scale as mp, it seems plausible that type la su- 
pernovae would occur within universes with significantly 
differing values of mp. This ensures the production of 
heavy elements within these universes. 

There is a possible caveat to this result. Within the 
context of another anthropic analysis, it has been re- 
marked [HI that the relatively low production of oxy- 
gen by type la supernovae may significantly hinder the 
formation of life if oxygen is not generated elsewhere. 
However, it is not clear that this suppression, roughly 3- 
8% relative type II supernova [6lj], is sufficient to render 
life overwhelmingly unlikely. The arguments of ref. fl8j 



were aimed against a scenario where type II supernova 
would definitely not occur. Since these supernova occur 
for a wide range of stellar masses within our universe, it 
is plausible that universes with mp not too unlike ours 
will also contain type II supernovae. It is beyond the 
scope of this paper to investigate more precisely for what 
values of mp these supernovae will occur. 



G. Stability of Stellar Systems 

As is illustrated in refs. fl5l . Il9j . an important an- 
thropic constraint derives from requiring that stellar sys- 
tems are stable against cosmic disruptions. Specifically, 
if a second star grazes too close to an existing stellar sys- 
tem, then a habitable planet may be thrown out of its 
anthropically fortuitous orbit. Here we seek a constraint 
to ensure that such encounters are typically too infre- 
quent to interfere with the evolution of life. First, we 
define a destructive encounter rate, 



7 ~ n+cr+Vp , 



(58) 



where n+ is the number density of stars, <t* = irb 2 is the 
cross section for an encounter with "fatal" impact param- 
eter b, and v p is the typical peculiar velocity of a star. 
Note that all of the stars within a given neighborhood 
have the same circular velocity; thus the circular velocity 
does not contribute to the encounter rate. 

The typical peculiar velocity of the stars in a galaxy 
is approximately determined by the temperature of the 
constitutive baryons during the phase of star formation. 
Since the baryons in a galaxy quickly cool to about Th ~ 
10 4 K and cool relatively slowly thereafter, we take this to 
be the relevant temperature. The corresponding peculiar 
velocity is then given by Eq. PS|) . For mp = 1 this 
gives a typical peculiar velocity of about v p ~ 20 km/s 
which agrees well with observation. The number density 
of stars n+ is equal to the number density of baryons 
divided by the typical number of baryons within a star, 
A*. Recall from Section Til El that the number density of 
baryons within a galactic halo is equal to fbPvir/mNfib- 
Therefore can be written, 



/*/fcPvir 

m N (j, b N+ ' 



(59) 



where /* is a fudge factor inserted to account for the 
increased density of the galactic disk, for the fraction 
of baryons that do not end up within stars, and for 
any clustering that may be involved in the star forma- 
tion process. A typical star within our universe contains 
10 57 baryons; therefore the results of Section III Fl sug- 
gest A* — 10 57 mp//t^ — 10 57 rhp/fib- To estimate /* is 
somewhat more challenging. 

The factor /+ accounts for several effects. For example, 
stars may form in clusters such that most stars exist in 
a neighborhood of higher density than the average den- 
sity of stars in a galaxy. On the other hand, a significant 
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fraction of baryons may compose a relatively diffuse in- 
terstellar gas and therefore not contribute to the stellar 
encounter rate. As it is beyond the scope of this paper 
to compute the mp dependence of these effects, we sim- 
ply treat them as being independent of mp. Meanwhile, 
we also expect /* to be proportional to the relative den- 
sity of the galactic disk to that of the baryons in the halo. 
Rather than concern ourselves with the specific geometry 
of the galactic disk, we study a simple model to obtain 
the mp dependence. We expect the factor /* to roughly 
scale like 



R :i 



-ffdisk-RLk 



(60) 



where Hd\ R y is the typical disk thickness. This can be 
written 15311 . 



(61) 



Meanwhile, to solve for i?disk we note that the circular 
velocity for stars is given by both by Eq. (|44[) and by, 



(62) 



Equating these expressions gives i?disk oc X 2 f b R. Fi- 
nally, putting all of this together gives, 



J* 01 \8 „.9 ' 



A 8 v 2 



(63) 



Note the strong dependence of /* on the spin param- 
eter A. The spin parameter is a stochastic variable with 
statistical properties related to those of £ oq . For ex- 
ample, the Milky Way appears to be characterized by 
A « 0.06 [g^l while typical galaxies ma y h ave A a fac- 
tor of two larger or smaller than this |54j |. This and 
other factors suggest that the factor /* may vary widely 
among galactic environments within any particular uni- 
verse. In addition, as explained above we have ignored 
several effects that might enter into /* . For concreteness 
we normalize /* to the value that describes our solar en- 
vironment in the Milky Way, /* ~ 10 5 [Hj]. This gives 



U ~ 10 5 mp" W b V /3 Pe q /3 CeqS(M) , 



(64) 



where the dependence of /* on mp stems entirely from 
its dependence on /& and T v ; r . 

It is left to calculate the impact parameter for fatal en- 
counters. We are specifically interested in the persistence 
of stellar systems that contain a planet in orbit about a 
star such as those considered in Section Hi Fl In addition, 
we focus on planetary orbits that receive electromagnetic 
radiation with an intensity that is comparable to that 
from the sun at the orbit of the earth. We then assume 
that an encounter will not be devastating to such a stellar 
system if the gravitational field from the grazing star is 
less than a tenth that of the primary star in the vicinity 
of the orbiting planet. Therefore we approximate b to be 



roughly v 10 times the radius of orbit for a planet receiv- 
ing about the same stellar intensity as the earth but in 
orbit about the stars studied in Section ["II Fl Note that 
the luminosity of a star is L oc R 2 Tf while the intensity 
at a distance r is I cx L/r 2 . Therefore 



10 r a 



RqT 2 i?g 



(65) 



where the subscript designates a quantity for the sun, 
<g> designates a star with surface temperature T® = 3500 
K in our universe, Rq is the radius of a star with this 
surface temperature within a universe with a different 
value for mp, and r au is one astronomical unit. Eq. (|65j) 
has been written so that every quantity can be evaluated 
within our universe except for Rq/R® rs mp, which is 
deduced using the results of Section III Ft 

We put all these results together to obtain 7. The con- 
straint that stellar systems typically survive a dangerous 
close encounter for long enough that life may evolve is 
Tevol 7 ; where again r GVO i = 5 x 10 9 yrs. In terms of 
cosmological parameters this is, 



3 ,-.-1/2 f-e d -4/3?-4 -2/3 - 
Jb Van Sen P A 



m P fi b 



4 > 10" 



(66) 



where we have used the models of ref. [59fl to substitute 



R®T^/R G T 2 r* ^Lq/Lq r* 0.14. Decreasing m P re- 
duces the cross section for dangerous impacts, since the 
'anthropically-favorable' radius decreases, yet increases 
the number density of stars. The net effect is a con- 
straint against decreasing mp. The explicit constraint 
implied for mp is given by 



g m p + g p 



2/3 



3+401+6/3,-2/3 



V 4 > 10~ 5 , (67) 



where fib depends relatively weakly on mp and has been 
ignored. The curves that saturate this inequality are dis- 
played in Fig. [T] using the label "7 _1 = r evo i-" 



H. Summary 

Let us now summarize the results of the previous sec- 
tions. Many anthropic constraints depend on the pri- 
mordial curvature perturbation £ eq and on the baryon to 
photon ratio rj. Lacking any standard model for the gen- 
eration of cither of these, we write them generically as 
£ cq w tfip" and r\ w rhp , where mp is the ratio between 
the apparent Planck mass and the value obtained within 
our universe. For the most popular models of inflation, 
a ranges between one (m^tp 2 chaotic inflation) and three 
(hybrid and natural inflation). Meanwhile, popular mod- 
els of baryogenesis give between zero (efficient lepto- 
genesis and electroweak baryogenesis) and 3/2 (specific 
models of SUSY Affleck-Dine baryogenesis). 

Most of the anthropic constraints under consideration 
are displayed in Fig.[TJ for representative values of a and 
j3. We assume a WIMP candidate to dominate the dark 
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FIG. 1: Anthropic constraints on rhp, plot as a function of halo mass scale /i, for a = 1, f3 = (top left panel); a = 1, j3 = 3/2 
(top right panel); a — 3, /3 — (bottom left panel); a — 3, (3 — 3/2 (bottom right panel). The logarithms are base ten and 
the empty circle corresponds to the mass scale of our galaxy with mp as observed within our universe. The region within the 
li — rhp plane that is excluded by any constraint is the region that does not include this circle. The parameters a and f5 as 
well as the labels on the curves are defined in the text (see for example the summary of Section [II Hp . Note that the galactic 
cooling constraint does not include the effects of molecular cooling. 



matter density. Many constraints depend on the total 
mass within the galactic halo for which they are evalu- 
ated. This mass is denoted /i and is measured in units 
of the Milky Way mass, or 10 12 solar masses. Note the 
empty circle in each panel of Fig. [TJ This corresponds to 
a mass scale equal to the mass of the Milky Way with a 
Planck mass equal to the value obtained within our uni- 
verse. The region within the /i-mp plane that is excluded 
by any constraint is the region that does not include 
this circle. For clarity we do not display the constraints 
that primordial black holes form a sub-dominant contri- 
bution to the energy density of the universe and that 
the dark matter dominates over baryonic matter. These 
are weaker constraints than those displayed in Fig. [1] and 
they are easy to calculate from Eq. (|20|) and Eqs. ([32]) . 
Finally, we note that many of the constraints in Fig. [1] 
are deduced by assuming that other constraints are satis- 
fied. For example, the "disk inst." curve is changed when 
the constraint represented by the curve "i rec = e~ 2 / 3 i V i r " 
is not satisfied. The continuous curves in Fig. Q] are in- 



tended to guide the eye. 

We annotate Fig. Q] as follows. A number of con- 
straints come from the various levels of structure forma- 
tion. The curve labeled "t coo i = r grav " (this curve has 
a distinctive "dorsal fin" shape) marks the separation 
between the mass scales of halos that contain baryons 
which cool faster than they (would) collapse and those 
that do not. As explained in Section Til E 21 this is one 
among a set of sufficient, but perhaps not necessary, con- 
ditions that allow for galaxy formation. Another one of 
these conditions is that galactic disks be Jeans-unstable, 
which occurs below the curve labeled "disk inst." in 
Fig. [T] Meanwhile, structure formation requires that 
over-densities separate from the cosmic expansion before 
the domination of the cosmological constant halts their 
growth. This requirement is filled below the curve la- 
beled "coo = 1.69." Finally, our analysis of structure 
formation assumes that galactic halos virialize after re- 
combination, which occurs for mp values located above 
the curve labeled "t rec = e~ 2 / 3 i v i r ." Alternative paths to 



16 



structure formation are discussed in Appendix |A"1 

We also consider a few anthropic criteria that are not 
directly related to structure formation. For example, an- 
other constraint that we consider is that galaxies not be 
so small that they are blown apart by internal super- 
novae. This will not happen if the binding energy of 
a galactic halo well exceeds the energy released via su- 
pernovae. This condition is satisfied for /i and mp to 
the right of the curve labeled "-Ebind = E SD ." In ad- 
dition, one might require that collisions between stellar 
systems be such that impact parameters so small as to 
dislodge a habitable planet occur on a timescale that is 
much larger than the evolutionary timescale, here taken 
to be Tcvoi ~ 5x 10 9 yrs. This constraint is satisfied above 
the curve labeled "7" 1 = r ovo i." Finally, one might wish 
to restrict attention to universes that contain stars that 
have surface temperatures in excess of about 3500 K and 
that have main-sequence lifetimes in excess of about four 
billion years. These correspond to positions above the 
line labeled "r* = r ovol " in Fig. [1] 

Except for the stellar lifetime constraint, every con- 
straint displayed in Fig. [T] depends on the size of the ini- 
tial over-density that eventually grows into a galaxy. The 
curves in Fig. [1] correspond to choosing this initial fluctu- 
ation to be the rms of the density perturbations at a scale 
/i evaluated at matter-radiation equality. However, the 
initial over-density describing any galaxy is a stochastic 
variable that may be larger or smaller than this. There- 
fore all of the curves in Fig. [1] will be shifted when one 
considers galaxies that are away from the norm. In ad- 
dition, the disk instability and close encounters curves 
(labeled "disk inst." and = r cvo i" respectively) de- 
pends very strongly on other stochastic quantities, such 
as the galactic spin parameter (see Section Hi G|) . There- 
fore the range of mp that is consistent with the above 
constraints is larger than the windows in Fig. [T] would 
suggest if one allows observers to arise in atypical envi- 
ronments within any given universe. 



III. THE PROBABILITY DISTRIBUTION FOR 
THE SCALE OF GRAVITY 

If the Planck mass mp scans across a landscape of uni- 
verses, then the value within any particular universe may 
not be uniquely determined. However, with an under- 
standing of the landscape and a calculus to regulate over 
a conceivably infinite number of infinitely expansive uni- 
verses, we may in principle calculate the distribution of 
mp. Since we cannot access any of the other universes 
within the landscape, such a distribution cannot be di- 
rectly tested. Nevertheless, we may still use this distri- 
bution to calculate the likelihood that we should observe 
the value of mp that we do. As we are forced to test this 
distribution using only our universe, we must be careful 
to account for any selection effects that would attenuate 
the distribution of mp . 

These selection effects generate a factor S that multi- 



plies the "prior" distribution 2. Thus we write the prob- 
ability to measure the Planck mass to be mp, 

P(m P ) = S(m P )I(m P ) . (68) 

The factor I(mp) may be taken as the likelihood for uni- 
verses with Planck mass mp to arise within the multi- 
verse, while S(mp) may be understood as the likelihood 
for observers to arise within those universes. We restrict 
our prior X to account for only universes exactly like ours 
except for their value of mp. This is equivalent to re- 
stricting the selection criteria in S. As mentioned in the 
introduction, there are many subtle issues that compli- 
cate the calculation of S and T. Our purpose here is not 
to resolve any of these issues based on technical grounds. 
Instead we explore an empirical constraint that may com- 
plicate some proposals to address them. 

One might expect the likelihood for a universe to sup- 
port observers to be proportional to the total baryonic 
mass within galaxies in that universe. (Note that here 
and below we do not presume proportionality factors to 
be independent of mp.) Meanwhile, the baryonic mass 
within galaxies is proportional to the total energy within 
a universe. This quantity diverges in proportion to the 
volume of the universe. Nevertheless, we may hope for 
a regularization scheme that allows for the volumes of 
universes to be compared. Since the energy density does 
not redshift during inflation, it is possible that when vol- 
umes are properly regulated, the ratio between the total 
energy densities of two universes will be proportional to 
the ratio of their inflationary expansion factors. 

While this argument may be intuitively appealing, 
such a prescription for volume-based weighting presents 
well-known difficulties [24| . For instance, its conclusion 
is crucially dependent on a specific global space-like slic- 
ing, which is ambiguous outside the horizon of any one 
observer. For example, an observer can chose a space-like 
slicing that is engineered to create a very large initial vol- 
ume for the observer's own universe, while also suppress- 
ing the initial volume of the universes of casually discon- 
nected observers. The suppression of an initial volume 
can be used to cancel the inflationary expansion factor, 
such that this slicing would give a dramatically different 
counting than the weighting described above. 

Nevertheless, in at least one proposal this ambiguity 
has been overcome and the result includes a selection ef- 
fect that weights universes according to their inflationary 
volume (see Garriga et al. (2006) among refs. [HI)- We 
assume that this result holds and write, 

S{m P ) = ^(mp)V(mp) . (69) 

Here A is proportional to the anthropic factor, which 
ultimately gives the likelihood per unit volume for some 
class of observer to arise within a universe. Depending on 
one's notion of an observer, A might include, for example, 
the baryonic mass fraction within galaxies, the fraction 
of stars with lifetimes in excess of few billion years, etc. 
The factor V is the inflationary expansion factor for the 
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universe 1 . 



V(m P ) 



3JV(m P ) 



(70) 



where N(mp) is the number of e-folds if inflation that 
typically occurs after a universe with Planck mass mp 
has arisen within the multiverse. 

If inflation is driven by a single canonical scalar field, 
then the number of e-folds of inflation is 



N 



1 



V 

rf *f> 



dip . 



(71) 



Here ff is the value of the inflaton when inflation ends, 
set by when the first slow-roll parameter equals unity, 
and fi is the value of the inflaton when inflation begins. 
Note that in general N depends explicitly on mp. As a 
specific example, consider chaotic inflation with an infla- 
ton potential V(f) = |m^</? 2 . This gives 



N 



4 mp 



(72) 



It seems evident that N will generically depend on mp. 
Yet without understanding the mechanism by which a 
universe is obtained within the multiverse, it is not clear 
what is the (typical) value of (fi and what is its depen- 
dence on mp. 

To illustrate that the explicit and implicit dependences 
of N on mp are not expected to cancel, three models to 
determine fi are now considered. The first model sets 
fi to be the value where classical evolution of f begins 
to dominate over the quantum fluctuations experienced 
whenever a mode exits the Hubble radius. This is set by 
the solution to 



1 y3/2 
12 7T m^V v 



1. 

If this is the case, the total number of e-folds is 



(73) 



(74) 



The second model assumes that fi is determined by 
where the inflaton energy density equals the Planck en- 
ergy density. In this case, 



N 



rr> 



2 ' 



(75) 



On the other hand, if fi is determined by where the 
inflaton energy density equals M 4 , then 



N: 



M 4 



m^mp 



(76) 



Not only does N generally depend on mp, but the 
dependence is very strong for mp near the value ob- 
tained within our universe. Consider for example the 
case of chaotic inflation with fi set by Eq. (|75j) . Then 



N ~ 10 5 mp, where mp is the Planck mass in units of 
the value obtained within our universe. Meanwhile, if 
Eq. (75|) sets the value of f. t , then N ~ 10 10 mp. Clearly 
different choices for fi , and in particular different models 
of inflation, will in general give a different dependence of 
N on mp. However, the dependence is always strong. 
This is because our universe experienced a large number 
(at least about sixty) of e-folds of inflation. 

The ambiguity over the mp dependence of N is not of 
concern. The important result is that so long as the de- 
pendence on mp of A and 2 is significantly weaker than 
the strong exponential dependence in V, then we expect 
mp to be most probably observed very near one of the 
boundaries of the anthropic range. The analysis of Sec- 
tion |TT] reveals that this is not the case with at least the 
factor A and the value of mp observed within our uni- 
verse. We illustrate this with an explicit example in Ap- 
pendix [5] There we show that even in the contrived case 
where N ~ 60 mp, the volume factor V(mp) overwhelms 
what appears to be one of the tightest anthropic con- 
straints. This pushes the expectation value for mp well 
beyond what we estimated to be the anthropic boundary, 
while the value obtained within our universe sits far down 
the tail of the distribution. This is exactly analogous to 
the runaway "c-problem" and the "Q catastrophe" in- 
troduced in refs. (28[. We refer to our example as the 
"mp-problem." 

The cr-problem and Q catastrophe were motivated by 
the fact that in many models of inflation the total number 
of e-folds of inflation depends on the inflationary param- 
eters that also set the level of density perturbations (the 
authors of refs. [28| use the notations a ~ Q ~ Q. For ex- 
ample, in chaotic inflation with potential V{f) = |m^ 2 
one finds N ~ C 1 - Therefore if the inflationary parame- 
ters may scan over the landscape, by the same argument 
given above we expect £ to be pushed to one of its an- 
thropic boundaries, whereas in our universe it sits com- 
fortably near the middle of the anthropic window [l5| . It 
has been pointed out [63| that this argument is not com- 
pletely satisfactory, since by hypothesis universes with an 
enormous number of e-folds are preferred. In such uni- 
verses, C may plausibly depend on different parameters 
during a long stretch of early inflation than it does near 
the end of inflation, when scales important to the forma- 
tion of structure are generated. Moreover, we note that 
the curvature perturbation is related to the first slow-roll 
parameter, 
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Inflation of longer duration requires a smaller ei, yet 
for inflation to end at all requires that at some point 
ej evolve toward unity. Therefore £ oq may be signifi- 
cantly decreased from its value during most of inflation 
by the necessary condition that ej interpolate between 
some very small value and unity by the end of inflation. 

We emphasize that the mp-problem is not hampered 
by these issues. That is, unlike C, mp is a constant within 
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any given universe. 5 We also emphasize that if the model 
of inflation that describes our universe exhibits an mp- 
problem, then allowing more parameters to vary across 
the landscape cannot mitigate this problem. That is, 
although allowing more parameters to vary might dra- 
matically shift the expectation value for mp after the 
additional parameters have been marginalized, this can 
only happen if the overwhelming majority of universes 
near the new expectation value have values for the other 
parameters that are very different than ours. We would 
still be left with the challenge to explain why we find our- 
selves in a universe like ours, and not with these different 
parameter values. 

There are significant caveats to this result. First of 
all, it is not clear that the selection effects in S should 
actually factorize as in Eq. (|69|) . Since the diverging vol- 
umes of sub-universes is one of the circumstances that 
complicates making landscape predictions, we cannot be 
assured that the resolution of this problem will result in 
universes with greater inflationary expansion factors be- 
ing more likely to harbor observers. Another caveat to 
this discussion is that little is known about the distribu- 
tion I(mp). As we have defined it, this term receives two 
separate contributions. One contribution comes from the 
distribution of mp values over the landscape, that is the 
frequency of mp values among the number of meta-stable 
states that are allowed by the underlying theory. A sec- 
ond contribution comes from the dynamics of the multi- 
verse, which may prefer certain meta-stable states over 
others as the multiverse evolves in time. This is because 
the tunneling and diffusion rates of quantum fields will in 
general depend on mp , such that meta-stable states with 
certain values of mp will appear more frequently within 
the multiverse than others. This mp dependence within 
I(mp) could be very strong; see for example the studies 
of quantum diffusion in refs. (23j . 

Therefore I(mp) may depend more sharply on mp 
than does V(mp), with a local peak within the anthropic 
range. This might at first seem incredibly fortuitous. 
However, the situation is very different from the case of 
the cosmological constant p\ . In that case we observe p\ 
to be very far from its 'natural' value and therefore we 
must presume a very diverse and densely packed land- 
scape in order for the value that we observe to exist at 
all. However, since we do not know the natural value of 
mp, its landscape window could be much smaller. 6 In ad- 



dition, it is possible that the landscape is not as densely 
populated as we have presumed, in particular once we 
restrict attention to meta-stable states in every way like 
ours except in the value of mp. For example, if the spac- 
ings between allowed values of mp are significant next to 
the size of the anthropic window, then our value of mp 
might be consistent with the shape of I(mp). 

Furthermore, the model of inflation that describes our 
universe may not actually exhibit an mp-problem. This 
would happen, for example, if the number of e-folds of 
inflation that describe this model were independent of 
mp or had a maximum for some finite value of mp. An 
interesting example of the latter case occurs when the 
effective Planck mass is not fixed within our meta-stable 
state, but evolves as in Brans-Dicke theory. This sce- 
nario has been studied in refs. 23], where it is shown 
that for some non-minimally coupled models of inflation, 
the inflationary expansion factor is maximized when in- 
flation ends at some finite value of mp. Note however 
that for this or any other model of inflation to avoid the 
mp-problem, it would have to generate more e-folds of 
inflation than all of the other anthropically viable possi- 
bilities within the landscape. Moreover, the value of mp 
that maximizes N would have to lie within the anthropic 
window. 

Finally, it is possible that the analysis of Section |TT] 
missed or underestimated an important anthropic condi- 
tion. This might appear as the most attractive possibil- 
ity, but one must be careful to appreciate the strength 
of the exponential dependence within V(mp). In order 
to cancel this exponential dependence and thus make the 
observed value of mp reasonably likely, an anthropic con- 
straint must appear to exponentially suppress the like- 
lihood for observers to arise within our universe. The 
observed prevalence of galaxies, long-lived stars, super- 
novae, and planets, along with the observation that our 
solar system does not seem to occupy a particularly over- 
or under-dense region within the Milky Way, all seem to 
suggest that this is not the case. Since there do not 
yet exist experimentally confirmed theories for inflation, 
reheating, and baryogenesis, it is still possible that one 
of these processes presents an anthropic selection effect 
that provides this exponential suppression. This possibil- 
ity is explored relative to the reheating temperature and 
baryogenesis in the context of the c-problem in ref. . 
Since the reheating temperature in general also depends 
on mp, this analysis applies equally to our scenario. 



5 Of course, the inflationary landscape hypothesis presumes that 
the fields <j> described in the introduction, c.f. Eq. JTJ, will evolve 
within the multiverse. Indeed, this is how the landscape is pop- 
ulated. However, we assume that the vacua defined by the fields 
tf> are selected prior to period of inflation in which we take in- 
terest, during which mp is constant. Specifically, the factor X 
is assumed to account for any selection effects due to the field 
evolution prior to this period of inflation, and the terms A and 
V are defined to apply only after a particular meta-stable state, 
with a specific value of mp , has been selected. 

6 It is tantalizing that within the context of weighting universes 



by inflationary expansion factors, chaotic inflation with N set 
by either Eq. H74H or Eq. 175 > pushes mp to larger values, while 
the largeness of mp relative to other mass scales is well-noted in 
our universe. Let us assume a fundamental scale M ~ Mqtjt- 
Then in this case we expect mp ^> Mguti arm it i s possible that 
mp ~ 10 3 A/gut is simply the largest that the landscape allows. 
Furthermore, Eq. 1741 pushes the inflaton mass m v to smaller 
values, and perhaps m v is the smallest that the landscape allows. 
Thus we obtain the apparent hierarchy mp 3> A/qtjt m v . 
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IV. ANTHROPIC CONSTRAINTS ON A AND 
THE SCALE OF GRAVITY 

It is straightforward to extend the analysis of SectionlHl 
to the case where both mp and the cosmological constant 
Pa may (independently) scan over a landscape. The only 
constraint that is affected by this generalization is the 
requirement that over-densities separate from the Hub- 
ble flow before their growth is halted by the domination 
of the cosmological constant. The maximum amplitude 
reached by a linear rms over-density in this scenario is, 
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(78) 



On the other hand, an over-density has separated from 
the Hubble flow when a linear analysis gives a > 1 .69 . 
Therefore an rms fluctuation will eventually form a halo 
if tJoo > 1.69, which gives the generalization of Eq. 



(79) 



where we find it convenient to henceforth use a eq in- 
stead of CeqS. This is the only result from Section |TT] 
that changes when p\ may scan over the landscape. 

Clearly, Eq. (|T9"|) is weakened as p\ is decreased from 
the value it obtains within our universe. In this case, 
Eq. ([28]) eventually ceases to be the strongest constraint 
and mp is bounded from above by one of the other curves 
in Fig. [T] We may also interpret Eq. (|T9")) as an upper 
bound on pa for a specified value of mp. In universes 
with a larger value of mp , p\ is then more tightly bound 
than in our universe. However, in universes where mp 
is smaller than in our universe, the bound on p\ may 
be significantly weakened. This effect can be dramatic. 
For example, if we take mp = 0.1 and if a = 1 and 
(3 = 0, then p\ may be increased by a roughly a factor 
of ten million and still satisfy Eq. (|T9|) . Of course, to 
determine the most likely range within which to observe 
Pa requires to determine the prior distribution I(pa, mp) 
and to incorporate all of the selection effects into a factor 
6>(pA,mp), as described in Section ITTT1 Both of these 
tasks are beyond the scope of this paper. 

Nevertheless, it is worthwhile to proceed but within a 
very simplified picture. While our level of analysis does 
not permit even an approximate landscape prediction, 
our results do imply restrictions on the dependence of S 
and X on mp . Our first assumption is that the landscape 
is so densely packed that we can approximate the prior 
distribution I(mp, pa) to be a continuous and smooth 
function of both mp and pa- Then we can write the 
probability distribution for p\ in the form 



P(Pa) oc / 5(mp,pA)2"(m P ,pA)rfmp . 



(80) 



We discuss in Section IIIII and in Appendix [B] how our 
universe appears extremely unlikely to be observed if S 



contains a factor proportional to the inflationary expan- 
sion factor. Since we wish to expose additional restric- 
tions on S and I, we now assume that S does not contain 
this factor. 

It is helpful to first consider the distribution P(pa) 
when mp is fixed to the value obtained within our uni- 
verse. This corresponds to taking S oc (5(mp — 1) and 
thus eliminating the integral in Eq. ([50)) . Refs. [7j argue 
that it is appropriate to restrict attention to only positive 
values of pa and to take the distribution X to be roughly 
independent of pa- Although the rms fluctuation cr cq is 
constrained by Eq. (|79p . any particular over-density may 
be larger or smaller than a cq . This implies that galaxies 
of a given mass will form in universes even when pa is 
larger than what is allowed by Eq. ([79"]) . On the other 
hand, galaxies of a given mass become statistically rarer 
as pa is increased. To account for this, it is customary 
to speculate that the likelihood for a particular universe 
to be observed is proportional to the fraction of its to- 
tal mass that collapses into galaxies with masses above 
some minimum p m i n 0]. This minimum galaxy mass is 
presumably set by other anthropic considerations. 

The spectrum of density perturbations is at least ap- 
proximately described by Gaussian statistics. Therefore 
a randomly selected co-moving volume may or may not 
collapse, depending on the size of the matter over-density 
contained within the volume. We parameterize volumes 
using the mass p that they enclose, measured in units 
of the Milky Mass, 1O 12 M . Then the likelihood that a 
mass p will eventually separate from the cosmic expan- 
sion is given by the Press-Schechter function [49] , 
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The percentage of over-densities that eventually virial- 
ize is a function of the enclosed mass p because the 
rms amplitude of the initial density perturbations <7 eq 
depends on p (see Section III E|) . The fraction of galax- 
ies that have mass between p and p + dp is (dF/dp)dp. 
Since F{p — > oo) = 0, this means that the fraction 
of mass contained within galaxies with mass above the 
mass scale p is simply the Press-Schechter function F 
evaluated at p. When only the cosmological constant 
scans over the landscape, p eq = <7 e q = 1- This gives 
P{pa) oc S(pa) oc F(p min ,p A ) @, Ell- 
in order to study the scenario where both pa and mp 
scan over the landscape, we adopt a very simplified pic- 
ture. First, we assume that X is independent of both pa 
and mp over the anthropically allowed window. We em- 
phasize that, unlike the case with pa, we are unaware of 
any physical justification for this assumption regarding 
mp. Second, we restrict our attention to galaxies with 
masses near the mass of the Milky Way. We perform this 
restriction simply so that we may ignore the scale depen- 
dence of anthropic constraints. It turns out that values 
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of rap somewhat larger than our own do not contribute 
significantly toward P(p\). To highlight this result we 
simply neglect all constraints on increasing mp. On the 
other hand, the selection effects that bound mp from be- 
low are very important when determining P{p\). For 
simplicity we consider selection effects from only one ad- 
ditional constraint; which is that stellar encounters are 
rare enough on average to allow for life to evolve in the 
intervening time. According to Fig. [T] this is usually the 
strongest constraint on decreasing mp. The exception 
appears to be the case of low a and low (3, where the 
galactic cooling constraint can interfere and the stellar 
lifetime constraint is not far below the close encounters 
constraint. We simply ignore the cooling constraint and 
note that we could just as well evaluate P{pa) for galaxy 
masses somewhat below the mass of the Milky Way to 
obtain a similar result. To account for stellar lifetimes, 
we impose a hard cut-off below mp = 0.1. 

As mentioned above, the mass fraction within galaxies 
with masses between p and p + dp is SF = (dF/dp)dp. 
This quantity in general depends on the time at which 
one looks at the universe. We count galaxies in the in- 
finite future, which is practically equivalent to counting 
galaxies at any time after the domination of pa- Then 
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The only p dependence within F stems from the de- 
pendence on Coo ex cr cq cx s(p), where s(p) is given by 
Eq. (|23|) . Within any given universe, to consider only 
galaxies with a particular mass p in the far future is 
equivalent to selecting only over-densities with a par- 
ticular amplitude at equality. This is because within 
that universe over-densities with smaller amplitudes will 
form galaxies with smaller mass while over-densities with 
larger amplitudes will form galaxies with larger mass (re- 
call that we look at the universe after p\ domination 
when the growth in over-densities has halted). The am- 
plitude of the initial over-density that is selected by look- 
ing at a particular galaxy mass p is the one that gives 
u(p) = 1.69 in the infinite future. 

We must now account for the close encounter con- 
straint mentioned above. This constraint is converted 
into a selection effect by noting that if the rate of dis- 
astrous encounters between stellar systems is 7, then 
the probability that a stellar system will survive for a 
time r is e _7T . The rate 7 is discussed in Section Til Gl 
Note that it depends on the amplitude of the initial over- 
density that seeded the galaxy. We restrict our attention 
to galaxies with masses near to the mass of the Milky 
Way. As described above, these galaxies only come from 
over-densities that satisfy a(p ps 1) = 1.69 in the infinite 
future. At equality, these over-densities have an ampli- 
tude 
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Now we can write the likelihood P ss that a stellar system 
will survive for at least a time r within this set of galaxies. 
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FIG. 2: The distribution P(pa) displayed against log(pA). 
The solid curve is P(pa) marginalized over universes in which 
mp may vary and with a — 1 and (3 = 0, the longer-dashed 
curve is the same quantity but for a = 3 and /3 = 3/2, while 
the shorter-dashed curve is P{ph) evaluated when mp is fixed 
to the value within our universe. All three distributions are 
for fixed galactic masses p = 1. See text for results obtained 
for a range of galactic masses. The normalizations are chosen 
for clarity. 



We take r 



5 x 10 9 yrs, which gives 
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where the dependence on a general mass scale p has been 
restored for future reference. 

So far our assumptions correspond to weighting uni- 
verses by the fraction of stellar systems that survive close 
encounters for longer than r evo i and that exist in galaxies 
with mass near to the Milky Way mass. We also require 
mp > 0.1 in order to ensure that sufficiently long-lived 
stars exist in these universes. Finally, we should account 
for the fact that the abundance of baryons relative dark 
matter will depend on the value of mp within each uni- 
verse. Putting all of this together gives the probability 
density, 



P(pa) oc / dp dm P f b (mp)P ss (m Pl p A ,p) 

J fimin Jo.l 

x-|-F(mp,pA,p) . (85) 

The full mp dependence of fb, P ss , and dF/dp is found 
by substitution of the results from Section [TTJ In general 
Eq. (|85|) is integrated over a window p nim < p < ^ max , 
but as motivated above, in our main analysis we restrict 
to a narrow window about p = 1. Finally, as explained 
previously, Eq. (|85|) makes the simplifying but unrealistic 
assumption that T(mp, pa) ~ constant. 

The results of a numerical computation of P(pa) are 
displayed in Fig. [2j For reference, we also display the 
result when mp is fixed to the value obtained within our 
universe (note that this corresponds to dF / dp\ ll= \ and 
not F(p = 1)). Our value of pa corresponds to the origin 
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on this graph. Evidently the assumptions of this section 
render the observation of pa at or below our value very 
unlikely. In fact, the fraction of P(pa) that sits below 
Pa = 1 is about 7 x 10~ 5 for a = 1, (3 = and about 
4 x 10 -4 for a = 3, j3 = 3/2. Since relatively large values 
of pa receive significant weight only when mp is relatively 
small, we also see that most of the weight of these distri- 
butions comes from values of mp that are below the value 
obtained within our universe. This is evidence of a sort 
of 'statistical pressure' that gives greater weight to those 
values of marginalized parameters that permit a larger 
value of pa- This is why it was unimportant to account 
for selection effects that constrain mp from above. 

Here we note the importance of the close encounters 
constraint and of recombination timing constraint (i.e. 
the "i roc = e~ 2 / 3 £ v ir" constraint) in bounding the an- 
thropically allowed variation of pa- Inspecting Fig. [1] 
indicates that of the constraints limiting mp from below, 
the close encounters bound is the strongest. For p, = 1, 
it bounds mp > 1CP 02 for a = 3 and f3 — 3/2, and 
approximately mp > 1CP ' 8 for a = l,/3 = 0. Inserting 
these limits on mp into Eq. (|79p gives maximum values 
for pa in good agreement with the peaks in Fig. Had 
only the star lifetime constraint mp > 0.1 been imposed, 
the maximum allowed value of pa would have been much 
larger. For a = 1 and (3 = 0, pA could be as large as 
4 x 10 7 . In the case of a = 3 and (3 = 3/2, pA can be as 
large as 2 x 10 14 . But the anthropic constraints in Fig. Q] 
for this latter model indicate that the "t roc = e~ 2 / 3 f v i r " 
constraint is stronger than the star lifetime constraint. 
Imposing the recombination constraint requires instead 
that pa be no larger than about 4 x 10 s . 

The analysis that leads to the curves in Fig. [5] gives 
at best a crude approximation for the actual probabil- 
ity distribution for pa ■ One improvement to the analysis 
would be to weight universes by the mass fraction that 
collapses into galaxies that have a range of anthropically 
favorable masses, instead of the fraction that collapses 
into only galaxies with the Milky Way mass. Includ- 
ing galaxies with greater masses will tend to push the 
weight of the distributions toward smaller pa, while in- 
cluding galaxies with smaller masses pushes the weight 
of the distributions toward larger pa- We have checked 
that under the assumptions outlined above, allowing for 
a range of galaxy masses 0.1 < fi < 10 tends to push the 
weight of the distributions P(pa) to slightly larger values 
of pa- 

Previous calculations of the distribution P(pa) inte- 
grate over all galaxy masses equal to or larger than the 
Milky Way mass. Although our anthropic considerations 
offer no reason to ignore galaxies with mass below that of 
the Milky Way, and galactic cooling constraints limit the 
formation of galaxies with larger masses, we nevertheless 
consider the evaluation of Eq. ([55)) for a range of masses 
1 < p < oo. Integrating over < pA < 1 gives the prob- 
ability P* that observers in such galaxies would observe 
a cosmological constant less than or equal to our own. 
Numerically, we find P* = 1 x 10~ 3 for a — 3,/5 = 3/2 



and P, = 2 x 10 -4 for a = l,j3 = 0. For comparison, we 
find that in our universe P* = 0.06. To reiterate, these 
calculations ignored any additional selection effects that 
might depend on p, such as the effects of different galactic 
cooling rates. 

Another improvement to the analysis would be to in- 
clude more mp-dependent selection effects. The abun- 
dance of heavy elements and of long-lived stars with 
appropriate surface temperatures both seem important 
when determining the likelihood for observers to arise 
within a universe. However, the analysis of Section |TT] 
does not shed light on how to calculate these selection 
effects. One thing that is clear is that the range of typi- 
cal stellar masses scales as mp. This means that galaxies 
of a fixed mass will contain more stars as mp is decreased. 
If all else were equal, as mp is decreased this would re- 
sult in a greater number of observers per unit baryon 
mass in a galaxy, which would tend to push the weight 
of P(/?a) toward larger values of pa- In addition, the 
rate of destructive encounters 7 is a function of stochas- 
tic variables, including for example the spin parameter A. 
The statistical distribution of these variables could tend 
to strengthen or weaken the close encounter constraint 
as a function of mp. However, when everything else is 
equal the 'statistical pressure' alluded to below Eq. ([55)1 
tends to give greater weight to those values of stochastic 
variables that allow for a smaller mp and larger pa- 7 

Of course, a proper calculation of P(pa) requires an 
understanding of the prior distribution for mp, X(mp). 
It must be emphasized that shape of I(mp) could dra- 
matically influence the shape of the distribution P(pa)- 
Therefore the results of this section are best understood 
as an empirical restriction on the dependence of I(mp) 
on mp. Since this is our main point in this section, let 
us be very explicit. The curves displayed in Fig. [3 sug- 
gest that within a very simplified landscape picture, it is 



7 We have confirmed this phenomena with the following simple 
example. The close encounter rate 7 is proportional to a factor 
/* that accounts for the increased density of the galactic disk 
relative the dark matter halo (see Section III E 3t . This factor 
depends sensitively on A, /* oc A -8 . In the preceding analy- 
sis, we normalized the factor /* so as to give the correct stellar 
density within our neighborhood of the Milky Way. However, 
A is a stochastic variable. N-body simulations suggest that the 
distribution for A can be approximated using 65j 

P\(X)d\ oc — exp [— 2 ln 2 (28.6A)] . (86) 
A 

This distribution has a peak at about A fs 0.03, while the Milky 
Way appears to be described by A 0.06 62]. This implies that 
typical values of A more tightly constrain mp than the value rep- 
resented in Fig. [2] Therefore, one might expect that when we 
treat A as a stochastic variable with distribution P\(\), that the 
weight of the distributions P(p\) will shift to smaller values of 
PA- in fact, the opposite trend occurs, as the previously men- 
tioned "statistical pressure" is such that the weight of the dis- 
tributions P(pa) actually shifts the location (in pa) of the peak 
by approximately an order of magnitude toward larger values of 
PA- 
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very unlikely to observe a value of p\ that is at or be- 
low the value within our universe. This means that if a 
landscape picture is to describe our universe, it should 
contain important ingredients that were neglected in our 
analysis. In addition, these additional ingredients should 
provide a strong emphasis for larger values of mp . Thus 
we conclude that for a landscape picture to describe our 
universe as among those that are likely to be observed, it 
is necessary that I(mp ) or some other neglected selection 
effect must receive the vast majority of its weight for val- 
ues of mp that are very near to or larger than the value 
obtained within our universe. We emphasize that the 
analysis of this section did not assume that the inflation- 
ary expansion factor enters into landscape calculations. 



V. CONCLUSIONS 

If the magnitude of the apparent Planck mass mp may 
scan across a landscape of possibilities, then there may 
exist universes with physical parameters and interactions 
in every way like those within our universe except for 
their value of mp. We have calculated the range over 
which mp may scan over such universes while still sat- 
isfying a number of anthropic constraints. Perhaps not 
surprisingly, if we combine all of the anthropic constraints 
we find a rather narrow window for allowed mp. The re- 
sults for WIMP dark matter and representative models 
of inflation and baryogenesis are displayed in Fig. [TJ Of 
course, the window for allowed mp is expanded if one 
loosens the anthropic criteria. 

More interestingly, this window will expand if an im- 
portant cosmological quantity is determined by a stochas- 
tic process. For example, many scenarios to generate a 
primordial curvature perturbation depend on the local 
vev of a light scalar field, as does the density of dark 
matter when it is determined by the axion. If these mod- 
els apply, then the curvature perturbation and/or dark 
matter density are not correlated with changes in mp, 
and a much larger window for mp may be able to satisfy 
anthropic constraints. Our purpose has been to calcu- 
late a minimal window for allowed mp, so we have not 
considered these possibilities in detail. 

Even a very small window for allowed values of mp has 
important implications for the landscape paradigm. In 
particular, the probability to observe a particular value 
of mp may be weighted by the inflationary expansion 
factor of universes that contain that value. This effect 
inputs a strong exponential dependence on mp into the 
probability distribution, which must be offset by another 
strong selection effect near the peak of the distribution. 
This other selection effect could be a very sharp peak or 
boundary to the underlying landscape distribution; oth- 
erwise the effect must come from an exponentially strong 
anthropic dependence on mp. Such a strong anthropic 
dependence on mp would be in conflict with the obser- 
vation that mp has even a narrow anthropic window in 
our universe. This is another example of the runaway 



inflation problem discussed in the recent literature. 

We also consider the anthropic window for the cosmo- 
logical constant A when both A and mp are allowed to 
independently scan over the landscape. Even when the 
allowed range for mp is relatively narrow, it still allows 
for a significant broadening of the allowed range for A. 
This is because A is only constrained by the necessity that 
cosmic structures separate from the Hubble flow before 
A domination. Meanwhile, the time at which structures 
separate from the cosmic expansion is proportional to a 
high power of mp. The result is that even for values of 
mp within the small allowed windows of Fig. [TJ A may 
be over ten million times larger in other universes than 
it is within ours. Just because A may be larger does not 
automatically imply that our value of A is less likely to 
be observed, since selection effects may ultimately weight 
smaller values of A more than larger values. We perform 
a very basic calculation which suggests that anthropic 
selection effects tend to make larger values of A more 
likely to be observed. This suggests that the observation 
of a cosmological constant at or below the level obtained 
within our universe is very unlikely unless unknown an- 
thropic selection effects or the underlying landscape dis- 
tribution of mp is dominated by values very near to or 
larger than the value obtained within our universe. 
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APPENDIX A: NON-STANDARD PATHS 
TOWARD STRUCTURE FORMATION 

One might wonder what are the constraints on struc- 
ture formation if we do not assume that dark matter dom- 
inates over baryonic matter, or that virialization occurs 
after recombination. If dark matter does not dominate 
over baryonic matter, then the evolution of over-densities 
in the dark matter does not significantly affect the evo- 
lution of over-densities in baryons. Without apprecia- 
ble dark matter potential wells, baryon over-densities do 
not grow (even logarithmically) until after recombina- 
tion. This is because in the era before recombination, 
the Jeans length for the tightly coupled baryons, 

Rj = ttH-'vs , (Al) 

where v s = l/v3 is the speed of sound prior to recombi- 
nation, is always larger than the Hubble radius. Growth 
therefore does not occur in either the radiation or baryon- 
dominated era until after recombination. Between recom- 
bination and the domination of cosmological constant the 
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evolution of over-densities may be approximated by 

a ps (a/a rec )a lec • (A2) 

The spectrum of fluctuations at recombination cr rec is 
scale dependent in the sense that it is constant for 
scales larger than the Hubble radius at recombination 
but rapidly decreases to zero as one looks at smaller dis- 
tance scales. This is because of the tight coupling be- 
tween baryon and radiation over-densities, and because 
the latter decay after they enter the Hubble radius. 

As in the standard picture, after the domination of 
cosmological constant over-densities will grow by a factor 
of 1.44 and then stop. Thus the maximum amplitude 
achieved by a linear analysis of an rms fluctuation is 

CToo ps 1.44 x (a A /a Tec )a rec ps 5 x 1CT 2 p^Ceq • (A3) 

Here we have used that cr roc ps 5 x 10~ 5 £ oq on scales 
larger than the Hubble radius at recombination, and that 
recombination occurs at a temperature T rec ps 3000 K, 
where we ignore the logarithmic dependence of T rec on 
top and r\. The formation of structure still requires that 
a linear analysis gives > 1.69 before the growth in 
over-densities is halted by the domination of cosmological 
constant. This gives the constraint, 

p rcc C4 > 5 x 10 5 . (A4) 

Eq. (IA4JI constrains top according to 

Qto/ + ^TOp 1 ^) to p 3q > 5 x 10 5 . (A5) 

This constraint is much stronger than the constraint it 
replaces, Eq. (|2"9"|) . 

Allowing for baryons to dominate the matter density 
of the universe may affect the other constraints in Fig. [T] 
in two ways. First, the halo density and background 
density at virialization, p v ; r and p*, are now reduced by a 
factor of 3 x 10~ 6 due to the difference between cr rec and 
cr C q. Second, structure formation only occurs on scales 
greater than the Hubble radius at recombination, since 
sub-horizon perturbations are suppressed. Ignoring the 
top dependence in T rcc , this implies a minimum halo mass 
set by the horizon mass at recombination, corresponding 
to a scale of roughly fx min ~ 10 6 . It can be shown that 
no value of top satisfies all of the constraints displayed in 
Fig. [1] after these effects have been included. Dropping 
the constraint that virialization precede recombination 
does not change this result. 

We now turn to the second assumption of Section lll El 
which is that recombination occurs at least an e-fold of 
expansion before virialization. To investigate what hap- 
pens when virialization occurs before recombination, we 
adopt the following simplified picture. Dark matter over- 
densities grow when they enter the Hubble radius, and 
we assume that they become non-linear and virialize as 
they would within our universe. However, growth in 



the baryon over-densities is hampered by their interac- 
tion with the photon Hubble flow before recombination. 
Therefore we approximate that baryons do not partici- 
pate at all in the over-densities of the dark matter and 
are rarefied relative the halo density as they follow the 
Hubble flow. 

Within this simplified model, the final baryon fraction 
within a halo will be at most about /h/187r 2 , and will 
decrease by a factor of e~ 3 ps 0.05 for each e-fold of ex- 
pansion between virialization and recombination. How- 
ever, it turns out that only two of the constraints that 
we consider depend significantly on the baryon fraction 
of the halo. These are the disk instability constraint of 
Section lll E 31 and the close encounters constraint of Sec- 
tion III Gl (The explicit /& dependence in the galactic 
cooling constraint of Section HIE 21 is canceled by an im- 
plicit dependence within A c .) To explore whether this 
situation opens a new window for allowed values of top , 
it is helpful to adopt the following picture. Instead of 
simply eliminating the recombination timing constraint 
of Eq. (|3Tj) . we continuously weaken it. For example, we 
may demand that recombination occur at most N Iec e- 
folds of expansion after virialization and then study the 
above constraints as N Tec is increased. 

When we do this, we find that the curves in Fig. [T] 
corresponding to the disk instability constraint, the close 
encounters constraint, and the recombination timing con- 
straint all slide downward as iV rec is increased. This shifts 
the allowed window for top such that larger values of 
top, including the value obtained within our universe, be- 
come excluded as lower values become allowed. It turns 
out that the disk instability curve slides downward at a 
faster rate than that of the recombination timing curve, 
so that as the window for allowed TOp moves to smaller 
top it also grows smaller. Ultimately, the window gets 
pushed against other constraints, such as the stellar life- 
time constraint or the galactic cooling constraint, and 
disappears. This happens at about N Tec ps a few. 



APPENDIX B: ANALYSIS OF A STRUCTURE 
FORMATION CONSTRAINT 

In Section IIIII it is argued that if the probability to 
observe a particular value of top is weighted in part by 
the inflationary expansion factor of universes that con- 
tain that value of top , then it is overwhelmingly preferred 
that top should be measured at one of the boundaries of 
its anthropic window. It is clear from the discussion of 
Section [TT] that our value of TOp is not at either of its 
anthropic boundaries. Nevertheless, it is worthwhile to 
investigate more quantitatively just how 'far' is our value 
of top from its anthropic boundaries. For simplicity we 
investigate the selection effect from only one anthropic 
constraint. Specifically, we look at the structure forma- 
tion requirement that halos virialize before the domina- 
tion of cosmological constant (Section HIE 1|) . Note that 
this provides the tightest constraint on TOp according to 
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1.2 



1.4 



1 . 6 




mp dependence is given by 



FIG. 3: The normalized distribution P(mp) for N — 60 mp 
(solid) and N = 601n(mp) (dashed) for when the landscape 
distribution for mp depends on the inflationary expansion 
factor. Our universe corresponds to trip = 1. In both cases F 
is defined using a = 3 and = 3/2. Although in both cases 
selection effects appear to prefer a specific range for mp, our 
value is far outside of this range. 



the curves in Fig. [T] 

The arguments of Section IIIII are appropriate primar- 
ily when the landscape is so densely packed that we can 
approximate the prior distribution I(rap) to be a contin- 
uous function of mp within the anthropic window. The 
probability to observe mp to lie within the range dmp 
can then be written, 



P(rrip)dmp oc A(mp)V(mp)2(mp)dmp 



(Bl) 



where the factors on the right-hand side are defined in 
Section lTTTl Of course, we are assuming that universes are 
weighted in part by their inflationary expansion factor V. 
Since we have no knowledge about the shape of I(mp), 
we take I(mp) » constant. As suggested above, we take 
the anthropic factor A to be conditioned by only the 
constraint that halos virialize before the domination of 
cosmological constant prevents this. 

To proceed, we assume that the likelihood to observe a 
given value of mp is proportional to the baryon fraction 
within galaxies of mass greater than or equal to the mass 
of the Milky Way, 10 12 M Q . Allowing for smaller galaxies 
or allowing for observers that do not require a galactic 
environment can only expand the window of allowed mp . 
The Press-Schechter function [4^] gives the fraction of 
matter that collapses into a galaxy of mass greater than 
or equal to a given scale. It is derived in Section [IVI and 
given by Eq. (|8l"j) . We reproduce it here for convenience, 



F = erfc 



0.373 



PA 



cq^cq 



1/3' 



(B2) 



The prefactor comes in part from evaluating F in the 
infinite future and at the Milky Way mass scale. The 



( PA 



VPcqO- 3 



cq 



1/3 



6 



6 



-4/3 



(B3) 



Note that when a and j3 are positive, F is a decreasing 
function of increasing mp . We are interested in the bary- 
onic matter within galaxies. Therefore A should contain 
a factor of the baryon fraction given by Eq. (fT5|) , along 
with F. 

Finally, we take V(mp) oc e 3N for N e-folds of infla- 
tion. We require that N(mp) be an increasing function 
of mp so that mp is pushed to larger values, saturating 
the constraint in Eq. (f2U)) . In addition, we want a = 3 
and = 3/2 so that this constraint on increasing mp 
is as strong as possible. Rather than propose a specific 
model of inflation, we assume that one can contrive a 
model with the relatively weak dependence N w 60 mp. 
Then putting all of our assumptions together gives the 
probability distribution, 



P(m P ) = Mf b (mp)F(mp)e 3N(m ^ , 



(B4) 



where Af is a normalization factor. The normalized dis- 
tribution P(mp) is displayed in Fig. [3] for a — 3 and 
(3 = 3/2. Although it is intriguing that in this scenario 
the erfc function overcomes the exponential volume fac- 
tor, this happens for a value of mp somewhat larger than 
the value that we observe. This 0(1) change in mp is sig- 
nificant due to the exponential sensitivity of P on mp. 
Indeed, the fraction of the distribution function P(mp) 
that sits below mp = 1 is completely negligible com- 
pared to that which sits above (explicitly, this fraction 
is roughly 10 -53 ). This distribution is so sharply peaked 
because N is relatively large. For example, expanding 
about the local maximum gives P ~ exp[— cA^A 2 ] where 
c ~ 0(1) and A is the difference between mp and its 
value at the maximum. Even if the number of e-folds de- 
pends very weakly on mp, for example N rs 601n(mp), 
we still find the preference for larger mp to be overwhelm- 
ing. 

We now provide a final point of clarification. A careful 
reader may notice that according to Fig. [31 the values of 
mp that are most likely to be observed lie well outside the 
anthropically allowed windows of Fig. [TJ This is because 
the relevant curves in Fig. [1] are calculated by assuming 
that all over-densities have initial amplitudes equal to 
the rms amplitude. Meanwhile, Fig.[3jtakes into account 
that the initial amplitude of an over-density is at least 
approximately a Gaussian random variable. The discrep- 
ancy between the results in Figs. [1] and [3] reflect that un- 
der the assumptions of this appendix, the overwhelming 
majority of galaxies stem from over-densities that begin 
with amplitudes many standard deviations away from the 
norm. Although these galaxies result from relatively un- 
likely initial over-densities, the fact that they arise within 
enormously larger universes more than compensates for 
this. This result stems from the sharp dependence on 
mp in the inflationary expansion factor. If the proper 
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landscape measure does not contain this factor, then the indicates, 
distribution for mp would be very different than Fig. [3] 
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